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ABSTRACT 
The  Center  for  Naval  Analyses  Computer  War  Game,  SEALIFT,  is  a  Monte- 
Carlo  simulation  designed  to  help  study  sealift  capabilities  in  an  ASW 
environment.   A  mathematical  model  of  the  SEALIFT  game  is  posited  to  ob- 
tain expected  value  results  approximating  those  of  the  SEALIFT  game.   The 
model  is  cast  in  the  Fortran  terminology  of  SEALIFT  and  an  effort  is  made 
to  accurately  reflect  the  SEALIFT  flow  chart  logic  in  its  development. 
Comparisons  with  SEALIFT  results  to  determine  the  model's  reliability  and 
accuracy  are  not  made. 
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1.   Introduction 

Objective.   In  1963-1964  the  Center  for  Naval  Analyses,  CNA,  of  the 
Franklin  Institute  constructed  a  rather  large  and  complex  computer  war  game, 
SEALIFT  I  1 1.  This  game  was  used  as  an  analytical  tool  for  the  Anti  Sub- 
marine Warfare  Study,  CYCLOPS  II.   Due  to  the  large  computer  running  time 
of  SEALIFT  an  after  the  fact  need  for  a  "quick  and  dirty"  approximation  to 
actual  SEALIFT  results  was  felt,  and  during  the  summer  of  1964  this  author 
was  privileged  to  work  with  the  CNA  CYCLOPS  study  group  on  that  approxima- 
tion.  This  thesis  is  an  outgrowth  of  that  summer's  work.   The  objective  of 
this  thesis  is  to  provide  the  logic  for  an  expected  value  analysis  of  that 
portion  of  the  5EALIFT  game  directed  toward  calculating  transport  ship 
losses  through  submarine  action.   This  thesis  is  then  a  "first  cut"  model 
to  approximate  a  model,  and  the  analysis  is  executed  under  simplifying 
assumptions. 

Background.   Functionally,  the  SEALIFT  game  pits  the  submarine  forces 
of  one  opponent  against  the  sea  transport  endeavors  of  another  opponent  in 
support  of  a  time  and  scope  limited  war  overseas.  The  major  game  output 
is  the  number  of  transport  ships  sunk  during  the  war  and,  through  inference, 
the  effect  of  these  losses  on  the  sealift  supply  rates. 

Structurally,  SEALIFT  is  an  event  store  Monte-Carlo  computer  simulation 
and  is  primarily  probabilistic  with  essentially  no  kinematics  or  geometry 
wired  in.   Although  certain  tactical  procedures  were  assumed  during  con- 
struction, they  may  be  circumvented  or  modified  by  the  input  parameter  set. 
This  set  is  on  the  order  of  four  hundred  and  contains  numerous  decision 
switches  as  well  as  numerical  data.   It  is  in  the  scope  of  this  parameter 
set  that  one  finds  the  flexibility  of  the  game  and,  unhappily,  the  bane  of 
forthright  analysis. 


The  focal  point  of  5EALIFT  is  the  submarine-convoy  encounter  and  of 
the  thirty-two  subroutines  which  comprise  5EALIFT,  five  form  the  substance 
of  this  event.   The  remaining  routines  are  directed  toward  making  the  game 
more  realistic  by  introducing  other  strategies  and  toward  controlling  the 
flow  of  the  game.   Figure  1  is  a  graphical  representation  of  the  game.   A 
general  and  more  detailed  description  of  the  game's  structure  and  logic  is 
quoted  from  5EALIFT  I   1  in  Appendix  I. 

Philosophy,,   All  submarine  loss  events  in  SEALIFT  except  those  incurred 
in  the  submarine-convoy  encounter  are  analyzed  here  utilizing  independent 
Bernoulli  trial  techniques  to  yield  simple  expected  value  results.   To 
effect  this  result,  several  event-time  dependent  routines  were  omitted  from 
consideration  and  a  few  were  modified  to  simplify  the  analysis,, 

Preliminary  investigation  indicated  that  the  submarine-convoy  inter- 
action was  the  most  sensitive  in  the  determination  of  transport  ship  losses 
and  that  simple  techniques  were  not  sufficiently  accurate  to  determine  these 
losses.   This  encounter  has  been  recast  in  the  framework  of  a  Markov  chain 
where  each  of  the  sequences  of  events  which  may  befall  a  submarine  upon 
convoy  passage  is  considered  a  state  of  the  Markov  chain.   Primary  effort 
has  been  directed  toward  accurately  effecting  the  Markov  structure,  and 
minimal  effort  toward  the  accurate  reflection  of  each  state  within  the  chain. 
Each  state,  i.e.,  sequence  of  events,  is  described  via  independent  Bernoulli 
trial  terminology,  again  with  approximations. 

This  analysis  has  been  made  with  some  trepidation.   It  may  be  argued 
that  detailed  Markov  analysis  is  pretentious  in  view  of  the  coarse  Bernoulli 
approximations;  nevertheless,  if  the  results  of  this  expected  value  analysis 
approximate  SEALIFT  results,  it  is  hoped  this  method  of  analysis  will  pro- 
vide new  insight  for  a  more  detailed  expected  value  analysis  of  the  entire 
SEALIFT  game. 


2.   Mathematical  Preliminaries 

Notational  usage  in  this  thesis.   The  following  symbols  are  used  in 

this  thesis: 

Symbol  Definition 

Pr(A)  The  probability  of  the  event  A. 

Pr(A/B)  The  probability  of  the  event  A  conditioned 

on  the  occurrence  of  the  event  Bo 

E(N)  The  expected  value  of  the  random  variable  N<, 

a  The  greatest  integer  less  than  or  equal  to  a< 

The  brackets  are  used,  additionally,  to  en- 
close a  matrix  when  its  elements  are  listed 
in  fullo 

(q. .)  The  matrix  whose  ij  th  element  is  q. .„ 

d 
A  =  B  A  is  defined  as  B,  or  A  is  identically  equal 

to  Be 
A  =  B  A  is  approximately  equal  to  B0 

n 

7(  A.  The  product,  A,  x  A_  x  <,.<>  x  A  „ 

.  ,   i  12         n 

i=l 

A.  The  sum,  A.  +  A_  ♦  . „ „  +  A  „ 

n 


r1  A.  I  he  sum,  A,  +  A„ 

L       1  12 

i=l 


iff  if  and  only  if 


Markov  chains.   It  is  generally  understood  that  if  a  physical  system 

can  exist  in  a  known  set  of  stable  states,  say  E..  ,  E_,  •••»  (finite  or 

countable),  and  if  with  each  state  there  is  associated  a  probability, 

Pr(E.)  =  p.,  that  the  system  is  in  that  state  (or  goes  to  that  state),  then 

the  probability  that  the  system  transits  the  states  E.,E.  ,...,E.   is 
n  J0   Jl       Jn 

the  product  IT      p .  •  The  basic  assumption  here  is  transitional  indepen- 
k=D   Jk 

dence,  i.e.,  the  chance  for  transition  from  a  given  state  to  another  does 

not  depend  upon  the  given  state  or  upon  past  states. 

Some  physical  systems  cannot  be  so  explained  but  require  the  dependence 

of  transitional  probabilities  upon  the  history  of  the  system.   The  simplest 

of  such  dependencies  is  that  the  next  state  in  sequence  depends  upon  the 

present  or  given  state  only.   In  this  case  the  outcome,  E.  ,  is  no  longer 

associated  with  a  fixed  probability  p,  ,  but  every  pair  (E.,  E.)  is  associ- 

k  i   J 

ated  with  a  conditional  probability,  p.  .,  which  denotes  the  chance  that 

the  system  goes  to  state  E.  given  that  it  was  in  state  E.  at  the  previous 

step.   In  general,  the  transitional  probability  p..  is  a  function  of  the 

sub-subscript  in  the  transit  E.  ,  ...,  E.  ,  and  the  process  is  called  homo- 

J0        Jn 

geneous  if  p . .  is  independent  of  this  sub-subscript.   (Only  the  homogeneous 

case  will  be  considered  in  the  expected  value  analysis  and  in  what  follows.) 

If,  in  addition,  the  probabilities,  p.,  of  the  initial  trial  i  are 

known,  the  probability  that  the  system  transits  the  states  E.  ,  E.  ,  ...,  E. 

J0   Jl       Jn 
n-1 
is  the  product  p .   x  77   P  •      ■>   The  conditional  probability  of  the 


J0   k=0   JkJk+l 


n-1 


transit,  given  the  initial  state  E .  ,  is  then   7T  P- 

J0  k=0  JkJk+l 

The  transitional  probabilities  can  be  arranged  in  a  matrix,  called  the 
transition  probability  matrix,  P 


(2.0) 


P  = 


Pll  P12  P13 


'21  P22  P23  '• 


P31  P32  P33 


where  the  standard  matrix  subscripts  ij  also  stand  for  transition  from  state 


E.  to  state  E..   Together  with  the  initial  distribution,  p.,  the  matrix 


completely  defines  a  Markov  chain  with  states  E   ,  E   ,  » <, .  *      It  should  be 
noted  that  the  elements  of  each  row  must  sum  to  one,  and  that  the  element 


p.  .  denotes  the  chance  for  transition  from  state  E.  to  state  E.  in  one  step<> 
ij  x  J 

The  natural  extension  is  to  investigate  the  concept  of  transition 
from  state  E.  to  state  E.  in  more  than  one  step,  say  in  n  steps „   This 
transition  can  occur  via  many  paths,  and  hence  the  conditional  probability, 
p.  .  ,  of  finding  the  system  in  state  E.  at  the  n+r  th  step,  given  that  it 
was  in  state  E.  at  the  r  th  step,  is  the  sum  of  the  conditional  probabilities 
for  each  possible  path.   In  particular, 


p!1'  .  p.. 


(2) 
ij 


■  LPaJP 


ir  rj 


The  equation 


(2.1) 


Pij 


(n)      r-     (n-1) 

-     Lp*.~ 


p 
xrrj 


and  the  Chapman-Kolmogorov  equation 


(2.2)     (ra+n) 
Pij 


V  (m)  (n) 
LP; 


p 
xr  rrj 


follow  by  mathematical  induction. 


William  Feller,  An_  Introduction  To  Probability  Theory  And  Its  Appli- 
cations  (New  York;   John  Wiley  and  Sons,  1957)  p.  348. 


The  probabilities,  p     ,  can  be  arranged  in  a  matrix  of  identical  structure 

as  that  given  by  equation  (2.0).   Coupled  with  equation  (2.1),  this  suggests 

matrix  multiplication  -  the  ij  th  element  of  the  n  step  matrix  is  the  inner 

product  of  the  i  th  row  of  the  one-step  matrix  with  the  j  th  column  of  the 

(n-1)  st  step  matrix,  i.e., 

Pn  =  P  x  P"-1 

Similarly,  equation  (2.2)  implies 
pm+n  =  pm  x  pn 

Also,  if  p .  is  the  probability  of  the  initial  state  E.,  the  unconditional 
probability  of  finding  the  system  in  state  E.  after  n  steps  is 

J 


(2.3) 


,(n) 


In)    <p 

j   =  ^  Pi 


x  p 


(n) 
ij 


A  particular  n  step  transition  matrix.   The  expected  value  analysis 
requires  the  n  step  matrix,  P  ,  formed  from  the  (one-step)  transition  prob- 
ability matrix,  P, 


(2.4) 


P  = 


Pl  P2  P3  P4  p5  P6  P7  P8 


Pl  P2  P3  P4  P5  P6  P7  PB 


00100000 

00010000 
00001000 

00000100 

00000010 

00000001 


=   (p..) 


where  p  =  p   =  p   of  the  previous  notation,  etc.   If  A  ,  A_,  A_,  A. 
are  defined  such  that 


Pl  P2 


P1P2 


,   A. 


P3  P4 


P3  P4 


,   A, 


p5  p6 


P5  P6 


.   Ay 


p7  PB 
P7P8 


then 


P  = 


and 
(2.5) 


Al  A2  A3  A4 
0   10  0 

0  0   10 

0  0  0   1 


and  in  general 


(AJL+I)A2    (A  +I)A3    (A  +I)A4 

0  0 

1  0 


0 

I 

0 

0 

0 

0 

An 
Al 

n-1 

(£     A*)A 
i=0     L 

0 

I 

0 

0 

0 

0 

n-1 

:£ 

i=0 


n— j. 


n-1 

<E  *J> 

i=<0 
0 

0 

I 


By  recursive  multiplication,  it  follows  that 


1  =   Pl  +  P2     X  1'   P  - 


and  therefore 


n-1 


n-1 

E 

i=l 


E  Ai  =   E  Ai  +  J  =   I!  <Pi  +  p2> 

i=0 


n-1 

E 

i=l 


i-1 


x  A   +  I 


1  ~  (Pjl  +  P2* 
1  "  (pl  +  P2J 


.n-1 


x  A  +  I 


Letting 


f2(n)   = 


1  -  (px  +  p. 
1  *  (Px  ♦  P2) 


,n-l 


;   fx(n) 


1  -  (px  +  P2) 
1  -  (Pl  +  P2) 


it  follows  that 


n-1 

E 

i=0 


Px  f2(n)  +1   p2  f2(n) 


Px  f2(n) 


p  f  (n)  +  1 


and  substitution  in  (2,4)  obtains 


(2.6) 


,        .n-1 
(Px  +  P2)    * 


Pl  P2 
Pl  P2 


6x2 


fx(n)   x 


P3  P4  P5  P6  P7  P8 
p3  p4  p5  p6  p7  P8 


'6x6 


3o   The  Model  for  the  Expected  Value  Analysis 

In  what  follows,  a  perusal  of  Appendix  I  and  of  SEALIFT  I   1   is 
assumed. 

Methodology.   The  heart  of  any  computer  game  is  the  flow  chart  logic, 
and  it  was  via  these  charts  that  SEALIFT  was  submitted  to  analysis.   Prior 
to  detailed  analysis,  the  game  was  examined  as  a  unit  to  cull  those  routines 
pertinent  to  the  submarine  from  those  of  a  bookkeeping  nature  and  other 
routines  of  non-interest,  including  those  which  determined  convoy  losses 
to  non-submarine  causes .   Of  the  thirty-two  routines  which  comprise  rthe 
SEALIFT  game,  twenty  were  discarded  as  irrelevant  to  the  objective  of  this 
thesis.   The  remaining  routines  were  again  scrutinized  for  syntax  as  related 
to  a  typical  play  of  the  game.   From  the  observed  functional  and  logical 
arrangement  of  these  routines,  the  all  important  abstraction  from  "real 
world,"  SEALIFT,  to  model  was  attempted. 

Concurrent  with,  but  not  independent  of,  this  abstraction,  a  variable 
which  would  be  definitive  to  the  expected  value  analysis  was  sought.   Bear- 
ing in  mind  that  the  motive  was  to  evaluate  losses  to  the  cargo  fleet  by 
submarine  offensive,  the  natural  and  obvious  choice  was  "the  total  expected 
number  of  torpedoes  fired  at  cargo  (vice  convoy)  ships  during  the  war."  As 
stated,  this  number  proved  to  be  unattainable,  the  primary  culprit  being 
the  time  and  event  dependencies  which  permeate  (are)  the  game.   It  appeared, 
however,  that  if  the  number  of  torpedoes  a  typical  submarine  would  fire  at 
cargo  ships  during  a  single  cycle  at  sea,  i.e.,  from  home  port  to  station 
to  home  port,  could  be  determined,  this  number  could  be  extended  into  a 
total  war  figure.   The  intuitive  appeal  of  this  single  cycle  approach  is 
that,  save  for  the  submarines  initially  on  station,  it  removed,  postponed 
really,  the  necessity  for  considering  the  element  of  time.   It  rendered 
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individual  submarines'  unique  and  time  dependent  lives  indistinguishable 
and  it  collapsed  the  immense  and  complicated  time-event  structure  into  a 
fairly  tractable  entity .   The  basic  cycle  is  represented  in  Figure  2.   The 
model  for  the  expected  value  analysis  developed  in  this  thesis  is  therefore 
the  logic  needed,  using  the  Markov  chain  process,  to  arrive  at  the  number 
of  torpedoes  a  submarine  would  fire  at  cargo  ships  during  a  single  cycle 
at  sea. 

In  developing  this  logic,  three  difficulties  became  immediately  apparent; 

1.  Aside  from  the  three  different  type  submarines  and  the  four  pos- 
sible sonar  boxes  to  which  they  could  be  assigned,  the  submarines  initially 
on  station  are  unique.   These  submarines  have  a  probability  distribution  on 
their  departure-from-on-station  time,  unlike  other  submarines  which  have  a 
fixed  mission  time.   This  probability  distribution  necessitated  a  modifica- 
tion of  the  associated  Markov  chain,  since  prior  SEALIFT  results  indicated 
this  cycle  to  be  a  very  sensitive  factor. 

2.  How  are  the  results  of  a  single  cycle  at  sea  "extended"  to  cover 
the  total  war?   In  particular,  what  assumptions  are  required  to  minimize 
the  effects  of  a  coarse  approximation?  How  are  excessively  small,  or  large, 
convoy  losses  to  be  interpreted? 

3.  How  can  the  event  and  time  dependencies  within  a  single  cycle  be 
handled  effectively? 

The  third  difficulty  was  resolved  by  essentially  ignoring  the  event 
and  time  dependencies  inherent  to  SEALIFT.   For  the  most  part,  these  depen- 
dencies proved  incompatible  with  either  the  cycle  abstraction  or  with  the 
Markov  chain,  or  both.   Where  possible,  however,  they  were  included  within 
the  framework  provided  by  this  thesis. 

The  second  difficulty  was  resolved  by  determining  the  minimum  number  of 
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cycles  a  live  submarine  would  make  during  the  war-   The  basic  cycle  was 
then  repeated  that  minimum  number  of  times  with  surviving  submarines  from 
one  cycle  serving  as  inputs  to  the  following  cycle.   Rationale  for  this 
choice  was  provided  by  two  facts.   First,  a  priori  information  from  the 
5EALIFT  game  indicated  that  where  the  submarine  was  almost  immediately 
overwhelmed,  the  significant  action  took  place  early  in  the  war.   In  this 
case,  the  number  of  cycles  would  prove  immaterial.   If,  however,  the  con- 
voys were  continually  decimated,  this  expected  value  analysis  would  show 
abnormally  high  convoy  losses,  and  this  would  actually  serve  the  purpose 
intended  by  5EALIFT. 

If  opposing  forces  were  fairly  balanced,  successive  cycles  would  show 
nearly  identical  convoy  losses,  again  providing  an  indication  of  the  progress 
of  the  war.   In  other  words,  under  the  assumptions  that  the  omitted  time  and 
event  dependent  events  are  not  significant,  the  analysis  provided  herein 
should  reflect  the  trends  of  a  corresponding  SEALIFT  game.   This  method 
of  extension  was  also  supported  by  the  fact  that,  almost  without  exception, 
the  omission  of  each  of  the  dependent  events  effects  an  over-estimation  of 
the  number  of  cargo  ship  losses.   The  choice  of  a  minimum  number  of  cycles 
would  then  partially  compensate  for  this  over-estimation. 

The  theoretical  resolution  of  the  first  difficulty  proved  trivial, 
falling  under  the  cognizance  of  the  integral  form  of  the  theorem  of  total 
probability.   The  evaluation  of  this  integral  proved  to  be  computationally 
messy.   The  functional  simplicity  of  the  probability  distribution  and  a 
linear  approximation  to  the  step  function  integrand  allowed  a  simple  but 
accurate  approximation  to  this  integral.   Each  of  these  difficulties  is  dis- 
cussed more  fully  in  subsequent  subsections,. 

Having  chosen  a  tentative  model  to  represent  the  game,  this  model  was 
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divided  into  its  component  parts  and  separately  analyzed.   The  basic  .cycle 
could  be  broken  into  two  phases,  (1)   the  submarine-barrier  penetration, 
and  (2)   the  submarine-convoy  interaction,  the  Markov  chain  of  this  develop- 
ments  The  probabilities  of  the  sequences  of  events  associated  with  each 
phase,  i.e.,  the  barrier  penetration  and  each  state  of  the  Markov  chain, 
were  analyzed  via  the  corresponding  5EALIFT  subroutines  or  parts  of  sub- 
routines.  The  SEALIFT  subroutine  logic  was  reformulated  in  the  structure 
of  a  tree  or  truth  table.   While  an  attempt  was  made  to  construct  the  trees 
to  match  the  subroutine  logic  exactly,  some  simplication  was  employed. 

Basic  assumptions  used  in  the  expected  value  analysis.   Perhaps  the 
most  definitive  assumption  made  in  the  expected  value  analysis  is  that  the 
submarine-convoy  encounter  should  be  the  significant  event  for  determining 
cargo  ship  losses.   By  implication,  then,  the  convoy  losses  to  air  attack, 
both  at  sea  and  in  port,  and  by  barrier  penetration  are  assumed  small  enough 
so  as  not  to  render  the  submarine  target-limited,  vice  torpedo-limited, 
and  also  that  each  convoy  be  initially  large  enough  to  comply  with  this 
restriction.   Thus,  Event  3  (Convoy  enters  barrier  C1C),  Event  4  (Air  attack 
against  convoy),  Event  6  (Convoy  enters  barrier  C2C),  Event  7  (Convoy  in 
delivery  port),  Event  8  (Convoy  enters  barrier  C2C  on  return  trip),  and 
Event  9  (Convoy  enters  barrier  C1C  on  return  trip)  were  classified  as  irre- 
levant to  the  analysis  for  this  reason. 

This  analysis  hinges,  then,  on  the  number  of  submarines  which  penetrate 
the  convoys,  and  prior  to  that,  the  number  of  submarines  which  reach  the 
subline.   For  this  reason,  all  the  submarine  loss  events  were  deemed  impor- 
tant.  In  the  subsequent  analysis,  however,  it  was  observed  that  the  train 
of  events  which  could  occur  to  one  submarine  was  dependent  upon  previous 
events  involving  other  submarines.   In  other  words,  the  stochastic  process 
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associated  with  a  particular  submarine  was  dependent  upon  the  realization 
of  the  process  associated  with  other  submarines  on  station  with  it  and  sub- 
marines which  proceeded  it  to  station.   Such  dependencies  are  inherently 
intractable  to  analysis  and  they  were  ignored  in  the  expected  value  models 
By  their  omission,  these  dependencies  introduced  a  cumulative  over-estimation 
of  the  number  of  cargo  ship  losses.   In  particular,  Event  17  (Submarine  enters 
variable  barrier)  calls  the  subroutine  VARDET  to  recalculate  the  reduced 
barrier-detects-submarine  probability  in  the  event  a  submarine  sinks  one 
or  more  of  the  barrier  ships.   Following  submarines  would  then  stand  a  better 
chance  of  slipping  through  that  barrier  until  the  loss  was  replaced  (Event  22 
(Replace  e  barrier  ship)).   A  concurrent  effect  is  the  torpedo  usage  on 
barrier  ships.   This  expenditure  directly  affects  the  ammunition  available 
to  be  used  on  cargo  ships  later  in  the  submarine's  cycle.   This  last  effect 
could  become  more  invidious  during  a  HUK  attack,  Event  15  (Submarine  versus 
HUK  group).   In  this  routine  the  only  ways  a  submarine  could  escape,  if  de- 
tected by  the  HUK  Group,  are  if  her  time  or  ammunition  runs  out  or  if  her 
torpedoes  miss  their  mark!   In  addition  to  introducing  event  dependencies, 
the  HUK  attack  event  presented  further  theoretical  difficulties  in  the  form 
of  time  dependencies.   Each  submarine,  while  at  sea,  is  subject  to  HUK  attack 
every  HUKHR  interval.   This  event  meshes  with  every  other  event  in  the  sub- 
marine's life  and  tends  to  hopelessly  confuse  expected  value  analysis. 

The  subroutines  FSEDET  and  ASWDET,  which  perform  the  same  function  for 
the  forward  screen  escort  and  the  ASW  screen  as  VARDET  does  for*  the  barrier 
ships,  are  called  during  the  play  of  Event  19  (Submarine  versus  convoy). 
Like  VARDET,  these  routines  influence  the  train  of  events  of  following  sub- 
marines.  Again,  parasitic  torpedo  usage  is  encountered;  however,  an  attempt 
was  made  to  approximate  this  usage  in  the  model  in  the  case  of  the  forward 
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screen  escort. 

Another,  and  time  related,  dependency  was  the  possible  monthly  change 
in  submarine-detects-convoy  probability  for  the  submarines  stationed  in  the 
fourth  sonar  box,  a  portion  of  Event  5  (Convoy  arrives  at  submarine  line). 
This  implied  an  incompatibility  with  the  cycle  concept.   Only  fortuitous 
combinations  of  the  number  of  cycles  of  the  war,  duration  of  the  war  and 
constancy  of  the  detection  probability  associated  with  the  fourth  sonar  box 
during  each  cycle  would  permit  easy  inclusion  of  this  dependency. 

Similarly,  Event  13  (Strike  at  submarine  source)  was  a  time  dependent 
event.   Provision  for  up  to  ten  air  strikes,  at  times  chosen  through  the 
input  parameter  set,  were  made  in  SEALIFT.   These  air  strikes,  too,  are 
incompatible  with  the  cycle  concept  used  in  the  model. 

To  summarize,  the  routines  ASWDET ,  F5DET ,  VARDET,  Event  15  (Submarine 
versus  HUK  group),  Event  22  (Replace  a  barrier  ship),  Event  13  (Strike  at 
submarine  source),  and  a  portion  of  Event  5  (Convoy  arrives  at  the  submarine 
line)  were  ignored  by  virtue  of  their  intractability  to  analysis. 

Finally,  Event  12  (Build  submarines)  provides  for  a  monthly  submarine 
building  rate.   While  this  event  is  also  time  dependent,  it  was  adapted  to 
the  model  by  applying  the  number  of  submarines  built  during  the  entire  war 
as  an  additional  input  to  the  first  cycle,  a  striking  illustration  of  the 
time  independence  of  the  cycle  concept. 

It  should  be  noted  that,  in  SEALIFT,  there  exist  at  least  two  levels  of 
logic,  the  internal  event  logic  which  computes  losses  to  convoys,  submarine 
losses,  etc.,  and  the  logic  which  controls  the  flow  of  the  game,  i.e.,  com- 
putes when  and  how  events  will  occur.   These  levels  are  frequently  meshed 
within  each  subprogram,  and  it  is  with  regard  to  the  first  named  level  that 
routines  were  ignored  or  classified  as  irrelevent.   The  second  type  logic 
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was  used,  implicitly,  whenever  necessary,  in  the  analysis.   In  addition, 
the  following  were  ignored  as  being  primarily  bookkeeping  type  routines: 
Event  1  (Build  ships),  Event  2  (Form  convoy),  Event  10  (Convoy  dispersed 
at  home),  Event  14  (Submarine  leaves  port),  Event  18  (Submarine  arrives  on 
station),  Event  20  (Submarine  goes  home),  Event  21  (Submarine  arrives  home), 
and  Event  24  (Replace  convoy  ships). 

The  events  from  SEALIFT  that  form  the  basis  for  the  analysis  are  there- 
fores   Event  5  (Convoy  arrives  at  submarine  line),  Event  11  (Initialize 
submarines),  Event  16  (Submarine  enters  constant  barrier),  Event  17  (Sub- 
marine enters  variable  barrier),  and  Event  19  (Submarine  versus  convoy K 

3.1  Markov  Chain  Structure  for  the  Expected  Value  Analysis 

For  a  given  submarine,  the  initial  Markov  chain  begins  after  that  sub- 
marines first  convoy  has  passed  the  subline.   The  initial  states  for  this 
Markov  chain  are  the  possible  sequences  of  events  which  could  have  occurred 
to  the  submarine  when  that  convoy  passed.   These  states  may  be  represented 
by  the  arborescence ,  with  root  "Submarine  on  Station,"  exhibited  in  Figure  3 
where  a  state  of  the  Markov  chain  is  defined  as  the  path  from  the  root  to  a 
terminal  point.   Letting  E.  denote  the  states  of  interest,  define 

E   =  submarine  does  not  detect  the  convoy. 

E_  =  submarine  detects  the  convoy,  does  not  penetrate  the  Forward 

Screen  Escort  (FSE),  and  is  not  killed. 
E   =  submarine  detects  the  convoy,  penetrates  the  FSE,  penetrates  the 

inner  screen  (ASW),  shoots  one  salvo  at  cargo  ships,  and  survives 

any  resulting  Search  and  Attack  Unit  (SAU)  action. 
E  =  submarine  detects  the  convoy,  penetrates  the  FSE,  penetrates  the 

ASW,  shoots  one  salvo  at  cargo  ships,  and  is  killed  in  resulting 

SAU  action. 
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SUBMARINE  15  ON  STATION 


SUBMARINE  DETECTS  CONVOY 


SUBMARINE  DOES  NOT 
PENETRATE  FSE 


SUBMARINE  IS 
KILLED  BY  FSE 


SUBMARINE  IS 
NOT  KILLED  BY 
FSE 


SUBMARINE  DOES  NOT 
DETECT  CONVOY 


SUBMARINE  PENETRATES  FSE 


SUBMARINE  PENETRATES 
ASW 


SUBMARINE  SHOOTS 
ONE  SALVO 


_L 


SUBMARINE 
IS  KILLED 
BY  SAU 


SUBMARINE  DOES 
NOT  PENETRATE 
ASW  (IS  KILLED) 


SUBMARINE  SHOOTS 
TWO  SALVOS 


SUBMARINE 
ESCAPES 


SUBMARINE 

IS  KILLED 

BY  SAU 


SUBMARINE 
ESCAPES 


FIGURE  3 
ARBORESCENCE  ASSOCIATED  WITH  THE  STATES  OF  THE  FIRST  MARKOV  CHAIN 
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E_  =  submarine  detects  the  convoy,  penetrates  the  FSE,  penetrates  the 
ASW,  shoots  two  salvos  at  cargo  ships ♦  and  survives  any  resulting 
5AU  action. 

E,.  =  submarine  detects  the  convoy,  penetrates  the  FSE,  penetrates  the 
ASW,  shoots  two  salvos  at  cargo  ships,  and  is  killed  in  the  re- 
sulting SAU  action. 

E_  =  submarine  detects  the  convoy,  penetrates  the  FSE,  and  is  killed 
by  the  ASW. 

Efl  =  submarine  detects  the  convoy  and  is  killed  by  the  FSE. 

This  Markov  chain  allows  for  the  possibility  that  if  a  single  salvo 
does  not  exhaust  the  torpedo  supply  of  a  submarine  during  an  attack,  there 
is  a  chance  that  the  submarine  will  be  unable  to  fire  again  at  that  convoy 
and  will,  therefore,  return  to  station  with  a  partial  load  of  torpedoes  to 
await  the  arrival  of  the  convoys.   It  is  assumed  in  this  model,  however, 
that  two  salvos  will  exhaust  a  submarine's  torpedo  supply.   This  assump- 
tion is  equivalent  to  the  assumption  that  the  necessary  train  of  events  which 
would  place  a  submarine  in  the  position  of  having  fired  more  than  two  single 
salvos  at  as  many  different  convoys  to  be  very  improbable,  and  neglects  this 
possibility.  i 

The  occurrence  of  states  E  and  E  imply  that  the  submarine  returns  to 
station  where  she  awaits  the  next  convoy  arrival,  which  is  to  say,  she  re- 
mains alive  in  this  Markov  chain.   State  E      is  followed  by  the  submarine's 
attempt  to  return  home,  and  is,  therefore,  an  absorbing  state.   Similarly, 

states  E  ,  E, ,  E.,,  and  EQ  are  absorbing  since  the  submarine  is  killed.   State 
4    fa    (       o 

E_  is  absorbing  for  another  reason;  it  is  the  entry  point  for  the  second 
Markov  chain  described  below. 

If  the  probability  of  the  initial  states  are  denoted  by  Pr(E.)  =  p., 
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the  one  step  transition  probability  matrix  is  given  by  equation  (2.4). 

If  a  total  of  n+1  convoys  pass  the  subline  during  a  submarine's  allotted 

on-station  time,  then  the  transition  matrix,  P,  must  be  raised  to  the  n  th 

power  since  the  initial  states  are  caused  by  passage  of  the  first  convoy  to 

yield  the  n  step  transition  matrix  P  .   The  resulting  matrix,  P  ,  is  given 

by  equation  (2.6).   This  matrix  and  the  initial  probability  distribution, 

p.,  suffice  to  determine  the  unconditional  probabilities,  P.   ,  that  a 
x  J 

submarine  will  be  in  each  of  the  possible  states  at  the  completion  of  n 
steps,  (n+1)  convoys  having  passed.   The  unconditional  probabilities,  P.   , 

J 

are  obtained  from  equation  (2.3)  and  from  them,  the  probabilities  of  certain 
relevant  events  can  be  obtained: 

(3.0)  P(A.)   =  Pr  (a  submarine  survives  this  Markov  chain  to  attempt 

passage  home/  n+1  convoys  have  passed  the  subline) 
.  P{n)  +  P*n)  +  P<n) 

(3.1)  P(A_)   =  Pr  (a  submarine  fires  one  salvo  at  the  cargo  ships  in 

this  Markov  chain/  n+1  convoys  have  passed  the 

subline) 

p(n)  .  p(n) 
3      4 

(3.2)  P(A_)   =  Pr  (a  submarine  fires  two  salvos  at  the  cargo  ships  in 

this  Markov  chain/  n+1  convoys  have  passed  the 

subline) 

p(n)  „  p(n) 
"5      6 

d 

(3.3)  P(B  )   =  Pr  (a  submarine  is  absorbed  in  state  E-  upon  passage 

of  the  r  th  convoy) 
=  p3  (p  +  p2)     ,   r  =  1,...,  n+1 
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The  second  Markov  chain  is  almost  identical  to  the  first  chain  and 
is  concerned  with  those  submarines  which  are  absorbed  in  state  E_.   These 
submarines  return  to  station  having  fired  a  single  salvo  of  torpedoes  at 
the  previous  convoy.   Assuming  that  a  total  of  (n+1)  convoy  cross  the  sub- 
line, the  number  of  steps  of  the  second  Markov  chain  depends  upon  which 
convoy  caused  a  submarine  to  be  absorbed  in  state  E_.   If  a  submarine  is 
absorbed  in  state  E_  at  the  passage  of  the  r  th  convoy,  the  (r+l)st  convoy 
"initializes"  the  second  Markov  chain  and  there  are  (n+1)  -  (r+2)  +  1  =  n-r 
steps  in  the  second  chain.   The  probability  that  a  submarine  enters  the 
second  Markov  chain  "just  after"  passage  of  the  r  th  convoy  is  given  by 
equation  (3.3). 

The  second  Markov  chain  is  described  by  its  states: 
F   =  submarine  does  not  detect  the  convoy,, 
F_  =  submarine  detects  the  convoy,  does  not  penetrate  the  F5E, 

and  is  not  killed. 
F-   =  submarine  detects  the  convoy,  penetrates  the  F5E ,  pene- 
trates the  A5W,  shoots  her  remaining  torpedoes  at  cargo 
ships,  and  survives  any  resulting  5AU  action. 
F   =  submarine  detects  the  convoy,  penetrates  the  FSE,  pene- 
trates the  ASW,  shoots  her  remaining  torpedoes  at  cargo 
ships,  and  is  killed  in  the  resulting  SAL)  action. 
F_   =  submarine  detects  the  convoy,  penetrates  the  FSE,  and  is 

killed  by  the  ASW. 
F^  =  submarine  detects  the  convoy,  and  is  killed  by  the  FSE. 

Denoting  the  probabilities  of  the  initial  states  of  this  Markov  chain  by 

Pr(F.)  =  q  ,  the  one  step  transition  probability  matrix  is  then 
i     i 
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Q  = 


qi  q2  q3  q4  q5  q6 


ql  q2  q3  q4  q5  q6 
0  0  10  0  0 

0  0  0  10  0 

0  0  0  0  10 

0  0  0  0  0  1 


=   (q..) 
ij 


and  the  n-r  step  transition  matrix,  by  analogy  with  equation  (2.6),  is 


,n-r 


,       «n-r-l 

(qj^  +  q2)     * 


ql  q2 


ql  q2 


i  -  (qx  +  q2) 
i  -  (qx  +  q2) 


n-r 


4x2 


q3  q4  q5  q6 


q3  q4  q5  q6 


4x4 


,  (n-r). 


If  R.   is  defined  as  the  unconditional  probability  that  the  submarine  is  in 
state  F.  after  n  steps,  then,  using  the  theorem  of  total  probability, 


(n) 


n+1 

I 

r=l 


Pr  (a  submarine  is  in  state  F.  after  n+1  convoys/ 
she  entered  the  second  chain  after  the  r  th 
convoy  has  passed)  x  Pr  (she  entered  the  second 
chain  after  the  r  th  convoy). 


n 

I 

r=l 


Q(.n"r)  x  P(Br)  ,   n  ^1 


where 


^-   r-*      (n-r)    ,       .r-1     > 

■   E    L  \  <ij     p3  (Pi +  p2)      •  n  = 

r=l  i=l 


(o)   d 

ij  Oiji    the   Kronecker  delta. 
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The  probabilities  of  the  remaining  events  of  relevance  can  now  be  ob- 
tained as : 

(3.4)       P(^a^      =  ^I    (a  submarine  survives  the  submarine-convoy  inter- 
action/ n+1  convoys  have  passed  the  subline) 


.  p(n>  +  p(„)  +  p(n)  ^   R(n)  +  R(„)  +  R(n)  +       +  ^n 


and 


(305)       P(Aj-)   =  Pr  (a  submarine  fires  her  remaining  salvo  in  the  second 

Markov  chain/  n+1  convoys  have  passed  the  subline) 

In  SEALIFT,  the  submarines  may  be  separated  into  two  sets,  those  which 
are  on  station  at  the  beginning  of  the  war  and  have  a  departure  time  distri- 
buted uniformly  on  an  interval,  and  those  submarines  which  arrive  on  station 
after  the  war  starts  and  have  a  fixed  mission  time.   Once  the  game  para- 
meters have  been  set  the  arrival  of  each  convoy  is  deterministiCo   As  a 
consequence,  the  variable  n+1  is  a  step  function  of  the  submarines  departure 
time  T,  i.e.,  n+1  =  g(T).   Representing  the  equations,  (3.0),  (3.1),  (3.2), 
(3.4)  and  (3.5)  symbolically  as 

Pr(A/n+l)   a  Pr  (of  the  event  A  given  n+1  convoys  have  passed) 
this  probability  becomes  Pr(A/g(T)  ).   Invoking  the  theorem  of  total  probabil- 
ity, Pr(A)  can  be  found  as 


■/ 


(3.6)       Pr(A)   =   /   Pr(A/g(t))  d  Pr(T^t) 

t 

This  equation  is  valid  for  both  sets  of  submarines  although  trivial  for 

those  with  a  fixed  mission  time.   For  the  submarine  with  a  fixed  mission  time 

Pr(A)  is  simply  the  integrand  of  this  equation  evaluated  at  the  upper  limit 
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of  t.   The  evaluation  of  equation  (3.6)  is  crucial  in  determining  the  ex- 
pected number  of  cargo  ships  sunk  per  cycle. 

This  completes  the  construction  of  the  Markov  chain  for  a  generalized 
submarine.   To  apply  this  structure  to  the  5EALIFT  game,  the  integrals  Pr(A) 
and  the  transition  probabilities  p. .  and  q. .  must  be  expressed  in  terms  of 
of  the  FORTRAN  variables  of  SEALIFT.   The  submarine  barrier  penetration  must 
be  analyzed  in  like  terms.   Finally,  the  evaluated  terms  must  be  synthesized 
to  answer  the  basic  question  "What  is  the  expected  number  of  torpedoes  fired 
at  cargo  ships  per  cycle?"  and  its  extension  "What  is  the  expected  number  of 
torpedoes  fired  at  cargo  ships  in  the  war?" 
3.2  Detailed  Analysis  in  Terms  of  the  FORTRAN  Variables  of  the  SEALIFT  Game. 

The  detailed  analysis  to  follow  will  be  executed  in  the  following  order, 
(1)  evaluation  of  the  integrals  Pr(A),  (2)  evaluation  of  the  barrier  penetra- 
tion, the  transition  probabilities,  and  related  expected  number,  (3)  syn- 
thesis to  determine  the  expected  number  of  torpedoes  fired  at  cargo  ships 
during  the  war.   The  following  list  of  5EALIFT  I  111  FORTRAN  variables  is 
pertinent  to  the  analysis,  and  the  following  convention  is  observed:   The 
letter  "I"  is  reserved  for  the  type  submarine  index;  however,  it  is  employed 
in  two  ways.   As  a  subscript,  e.g.,  ALTIM  (I),  "I"  can  take  on  the  values  1, 
2,  and  3.   As  a  part  of  a  variable  name,  e.g.,  NSUB'I'T,  "I"  is  set  off  in 
this  analysis  by  single  quotation  marks  and  can  take  on  the  values  1,  2  and 
N.   In  both  usages,  the  values  of  "I"  denote  conventional  number  1,  conven- 
tional number  2,  and  nuclear  submarines,  respectively. 

Variable  Name        Definition 

TBCON  Time  between  each  convoy  sailing  (hours). 

CONSPD  Speed  of  convoy  (knots). 
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CONHR 

ONEWAY 

ULTIME 
ALTIM(I) 

DSUBTA 

SUBSPD(I) 
NSUBTT 
FSUB'I'OS 
P'rBOX(J) 

NSUB'I' 

PBAR5C(I,J) 

PBARSV(I.J) 
PVBDET(I,J) 
VAIRDET(J) 


RGDT(K,I) 


SUBLINE 


TWAR 


Time  after  start  of  war  that  the  first  convoy  leaves 
home  port  in  hours. 

One  way  convoy  distance  from  home  port  to  delivery  port 
(n.  miles). 

Time  convoy  spends  in  delivery  port. 

Allowable  time  in  hours  for  type  I  submarine  before  it 
must  leave  its  station. 

Distance  from  submarine  home  port  to  submarine  station 
(n.  miles). 

Speed  of  type  I  submarine. 

Number  of  type  I  submarines  at  start  of  war. 

Fraction  of  type  I  submarine  on  station  at  start  of  war. 

Probability  that  a  type  I  submarine  on  station  will  be 
in  BOX  #  J,  J  =  1,  ..o,  4. 

Number  of  type  I  submarines  built  per  month. 

Probability  that  a  type  I  submarine  is  sunk  by  the  Jth 
constant  barrier,  J  =  1,  . ..,  4. 

Probability  that  a  type  I  submarine  is  sunk  by  the  Jth 
variable  barrier,  given  it  was  detected  by  this  barrier, 

w  —  X,   «..,   Ho 

Probability  that  a  type  I  submarine  is  detected  by  the 
Jth  variable  barrier,  with  barrier  at  its  full  quota  of 
ships e 

Part  of  detection  probability  of  variable  barrier  J 
that  is  not  supplied  by  ships  within  that  barrier  which 
might  be  sunk  by  a  submarine.   Thus,  if  a  mixed  barrier 
of  aircraft  and  ships  were  used  and  all  the  ships  were 
sunk,  the  probability  of  detection  drived  from  the  air 
barrier  (VAIRDET(J))  would  still  be  present.   The  "total" 
detection  probability  of  a  variable  barrier  is  PVBDET(I,J) 
+  VAIRDET(J).   For  a  ship  barrier  only,  VAIRDET(J)  should 
be  set  to  zero. 

Maximum  athwartship  detection  range  of  type  I  submarine 
against  convoy  type  K  sonar  conditions 

K=l  iff  good  sonar  conditions 

K=2  iff  poor  sonar  conditions 

Length  of  submarine  line  (n.  miles),  i.e.,  width  of 
patrol  area. 

Duration  of  war  (hours). 
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FAIRSC 

AT5FEU) 

PDFS(I,K) 


PDFER 


PDFCE(I.K) 
P5F5E(I) 
PFCE5(I,K) 
PSFCE(I,K) 

A5ALVQ ( I ) 
PDETU(I,K) 

VA(K) 


PBB(I,K) 
PKBB(I.K) 


Is  forward  escort  screen  made  up  entirely  of  aircraft 
(nan-sinkable)?   (YES=1,  N0=2h 

Number  of  times  a  type  I  submarine  can  try  to  repenetrate 
forward  screen. 

Probability  that  a  type  I  submarine  in  type  K  sonar  con- 
ditions (denoted  by  "(I,K)  submarine"  henceforward)  is 
detected  by  forward  screen  escort  while  penetrating 
screen  with  screen  at  its  full  quota  of  ships.. 

Part  of  detection  probability  of  forward  screen  that  is 
not  supplied  by  ships  within  that  screen  which  might  be 
sunk  by  a  submarine.   Thus,  if  a  mixed  screen  of  air- 
craft and  ships  were  used  and  all  the  ships  were  sunk, 
the  probability  of  detection  derived  from  the  air  screen 
(PDFER)  would  still  be  present.   The  "total"  detection 
probability  of  the  forward  screen  is  then  PDF5(I,K)  + 
PDFER.   For  a  ship's  screen  only,  PDFER  must  be  set  to 
zero. 

Probability  that  an  (I,K)  submarine  is  sunk  by  forward 
screen  escort,  given  submarine  has  been  detected  by  same. 

Probability  that  a  type  I  submarine  shoots  at  forward 
screen  escort  ship  in  an  attempt  to  break  contacto 

Probability  that  an  (I,K)  submarine  sinks  forward  screen 
escort  ship,  given  it  decides  to  shoot  at  it. 

Probability  that  an  (I,K)  submarine  is  sunk  by  forward 
screen  escort  ship,  given  this  submarine  has  already  sunk 
a  forward  screen  escort  ship,  i.e.,  flaming  datum. 

Number  of  anti-ASW  torpedoes  used  per  salvo  by  a  type  I 
submarine. 

Probability  that  an  (I,K)  submarine  is  detected  by  the 
ASW  inner  screen  with  direct  path  sonar,  given  it  was 
not  detected  previously  with  bottom  bounce  sonar» 

Conversion  factor  reflecting  improved  ASW  inner  screen 
detection  probability  with  direct  path  sonar,  in  type  K 
sonar  conditions,  because  of  previous  bottom  bounce  de- 
tection. 

Probability  that  an  (I,K)  submarine  is  detected  by  the 
ASW  inner  screen  with  bottom  bounce  sonar,, 

Probability  that  an  (I,K)  submarine  is  killed,  given 
that  it  is  attacked  or  re-attacked  following  a  detection 
with  bottom  bounce  sonar. 
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PAT(I,K) 
PCC(I,i) 


PABB(I,K) 


PKBBG(I) 


Probability  that  an  (I,K)  submarine  is  attacked  given 
that  it  has  been  detected  with  bottom  bounce  sonar  only, 

Probability  that  a  group  from  the  ASW  inner  screen  re- 
attacks  a  type  I  submarine  at  the  i  th  opportunity  in 
action  resulting  from  a  bottom  bounce  detection, 
i=l,  2,  30 

Probability  that  an  (I,K)  submarine  attacks  a  group 
from  the  ASW  inner  screen  in  action  resulting  from  a 
bottom  bounce  detection,, 

Probability  that  a  type  I  submarine  sinks  a  ship  in 
the  group  from  the  ASW  inner  screen  in  action  resulting 
from  a  bottom  bounce  detection. 


QTC  177 


PKINS(I) 


PSCA(I,i) 


PAA(I) 


PKSC(I) 


PSUBAT(I) 


PKOSC(I) 


Designator  for  possibility  that  an  engagement  between 
the  ASW  inner  screen  and  a  submarine  is  terminated  after 
an  unsuccessful  submarine  counter  attack,  e.g.,  due  to 
quiet  torpedoes  used  by  submarine  (QTC  equal  to  zero 
implies  engagement  is  to  continue;  QTC  equal  to  one 
implies  engagement  is  terminated.   No  other  values  are 
possible) . 

Probability  that  a  type  I  submarine  is  sunk  by  the  ASW 
inner  screen  during  action  following  a  detection  with 
direct  path  sonar. 

Probability  that  an  ASW  inner  screen  re-attacks  a 
type  I  submarine  at  the  i  th  opportunity  in  action  re- 
sulting from  a  detection  of  sub  by  ship  with  direct 
path  sonar,  i=l,  2,  3. 

Probability  that  a  type  I  submarine  attacks  the  ASW 
inner  screen  after  an  unsuccessful  attack  on  the  sub- 
marine by  the  screen  resulting  from  a  detection  with 
direct  path  sonar,, 

Probability  that  a  type  I  submarine  kills  an  ASW  inner 
screen  ship  in  action  resulting  from  detection  of  sub 
by  ship  with  direct  path  sonar. 

Probability  that  a  type  I  submarine  attacks  the  ASW 
inner  screen  following  an  unsuccessful  bottom  bounce 
detection  attempt  by  the  screen. 

Probability  that  a  type  I  submarine  sinks  an  ASW  inner 
screen  ship  without  prior  detection  of  the  submarine  by 
the  screen. 


QTU  398 


Designator  for  radiated  noise  level  of  enemy  anti-screen 
ship  torpedoes.   (QTU  equal  to  one  implies  a  quiet  tor- 
pedo; QTU  equal  to  zero  implies  a  noisy  torpedo.   No 
other  values  are  possible). 
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VM(K) 


CSALVO(I) 
PD5AU(I,K) 

PSSAU(I,K) 
ATTSAU(I) 
RASW(I) 
TORPS(I.ASW) 


Conversion  factor  reflecting  change  in  the  ASW  inner 
screen  detection  probability  with  direct  path  sonar  in 
type  K  sonar  conditions  because  of  flaming  datum  result- 
ing from  sinking  of  screen  ship  by  previously  undetected 
submarine 0 

Number  of  anti-shipping  torpedoes  used  by  type  I  sub- 
marine., 

Probability  that  an  (I,K)  submarine  is  detected  by  the 
5AV,  given  the  submarine  has  sunk  one  or  more  cargo  ships 
in  the  convoy  <> 

Probability  that  an  (I,K)  submarine  is  sunk  by  the  SAU, 
given  the  SAU  detected  the  submarine,, 

Does  a  type  I  submarine  attack  a  number  of  the  SAU  group, 
given  it  was  not  sunk  by  this  group,  YES=1. 

Number  of  times  a  type  I  submarine  can  re-attack  the 
convoy  after  SAU  action. 

Number  of  anti-ASW  torpedoes  per  type  I  submarine. 
Note:   If  only  one  type  of  torpedo  is  used  for  the  play 
of  the  game  and  it  is  desired  to  fire  this  torpedo  at 


all  targets,  then  it  should 
(I, ASW). 


must  be  stored  in  TORPS 


3.2.1  Evaluation  of  the  Integrals  Pr(A) 

To  evaluate  the  integrals  Pr't~A),  the  equivalences  between  the  first  and 

second  Markov  chains  are  exploited  to  condense  the  notation.   The  states  E  , 

E_,  E7,  and  E  are  seen  to  be  equivalent  to  the  states  F  ,  F  ,  F  ,  and  F, 
£    f       o  12   5       6 

respectively,  and  therefore,  p  =  q  ,  p  =  q  ,  p  =  q  and  pQ  =  p-„   For 

1    1   c  i.         (  3      ob 

s  =  (Px  +  P2) 


and 


1-s 


n+1 


f(n)   = 


1-s 


and  letting  D  denote  the  derivative  operator  with  respect  to  s,  the  evalua- 
tion of  the  unconditioned  probabilities,  P.   and  R.   yield  the  following 
relationships  s 
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o(n) 

pl  -  Pls 

D(n) 

P2  =  P2» 

P.  =  p  .f (n) ,  j  =3,  .00 ,  8 

Q(n)  n  n 

R1  =  Pl  p3  Ds 

R^n)  =  q2  p3  Df(n) 

R4P)  =  ^4  P3  Df(n) 

R^n)  =  p?  p3  Df(n) 

^^  =  P8P3Df(n) 

Using  these  values,  equations  (3.0,  (3.1),  (3.2),  (3.4)  and  (3.5)  from 

section  3.1  are  evaluated  to  yield 

(3.0)'   P(A,)  =  sn+1  +  p.  f(n) 


(3.1)   P(A2)   =   (p3  +  p4)  f(n) 


(3.2)   P(A3)   =   (p5  +  p6)  f(n) 

(3.4)'   P(A4)   =   (1  +  P3  D)  sn+1  +  (p5  +  p3q,3  D)  f(n) 


(3.5)   P(A5)   =   (q3  +  q4)  p3  D  f(n) 


Generically  each  of  the  integrals  Pr(A) 

Pr(A)   =    /.    Pr(A/g(t))   d  Pr(T  f  t) 


is  a  function  of  the  integral  N 

N   =  Jt        s9(t)  d  Pr(T^  t) 
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and  the  first  derivative  of  N  with  respect  to  s,  assuming  the  integration  and 

differentiation  operations  can  be  interchanged.   The  distribution  Pr(T  =  t) 

is  non-trivial  only  for  the  submarine  initially  on  station  and  is,  in  SEALIFT, 

uniform  on  the  interval  (0,  ALTIM)  ,  i.e«,   dPr(T  ^  t)  =  ..  _T,. 

ALI  1M 

To  explore  the  structure  of  the  function  g(t)  and  remain  within  the 

FORTRAN  concept  of  SEALIFT,  two  new  FORTRAN  variables  are  defined 

TCON(I)   =  The  game  at  which  the  I  th  outgoing  convoy  crosses 

the  subline, 
and 

RTCON(I)   =  The  game  time  at  which  I  th  returning  convoy  crosses 

the  subline. 
For  those  submarines  initially  on  station,  g(t)  is  zero  for  t  «=  TCON(l),  and 
steps  to  one  at  t  =  TC0N(1)„   It  then  increases  by  one  after  each  interval 
TBCON  until  the  first  returning  convoy  crosses  the  subline,  at  game  time 
RTCON(l),   Until  the  SEALIFT  game  parameters  have  been  set,  the  fine  structure 
of  g(t)  is  not  well  defined  after  RTCON(l).   It  does,  however,  increase  by 
two  every  TBCON  hours  after  RTCON(l).   Clearly  then,  the  integrals  Pr(A) 
involve  extensive  summations.   To  avoid  these  summations  a  piecewise  contin- 
uous linear  approximation,  h(t),  to  the  step  function,  g(t),  is  sought  such 
that  for  g(t)  =  h(t),  the  value  of  N  is  easily  found,  viz,, 


■/ 


,  +  .  h(t) 

(3.7)  N  =   /   sgU)   dt  I  •  -*- 


Ln(s)  dh(t) 
*  dt 


Incorporating  the  SEALIFT  FORTRAN  variables  and  the  fact  that  in  SEALIFT 

the  subline  is  positioned  half  way  between  convoy  home  port  and  delivery  port 

it  can  be  seen  that 
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TC0N(1»      "     2  "cONSPD     +  C™HH 

3  ONEWAY 

RTCON(l)      =    y     x     c^lpl      +     CONHR      +     ULTIME 

and   further, 

TCON(I)   =  TCON(l)  +  (1-1)  TBCON 
RTCON(I)   =  RTCON(l)  +  (1-1)  TBCON 
If  K  is  defined  as  the  number  of  outgoing  convoys  which  pass  the  subline 
before  the  first  returning  convoy  passes  the  subline,  then  K  is  the  largest 
integer  which  satisfies  the  inequality 

TCON(K)  %      RTCON(l) 


or 


(K  -  1)  TBCON  +  TC0N(1)  %     RTCON(l) 


K  = 


RTCON(l)  -  TC0N(1) 


TBCON 


+  1 


where  the  brackets  denote  "the  greatest  integer  less  than  or  equal  to." 
Introducing  a  new  FORTRAN  variable  H(I)  for  the  expected  value  analysis,  the 
values  of  H(I)  and  the  approximation  h(t)  are  defined  as: 
H(l)   =  h(0)   =  0 

H(2)  =  h(TC0N(l)  -  0)   =  0 

H(3)  =  h(TC0N(l)  +  0)   =  -J- 

H(4)  =  h(RTC0N(l)  -  0)   =  i+YBCON   (RTC0N(1)  -  TC0N(1) ) 

H(5)  £  h(RTC0N(l)  +0)   =  K  +  £ 

H(6)  =  h(ALTIM)   =   (K  +  |)  +  Tfi^0N   (ALTIM  -  RTCON(D) 

Figure  4  depicts  both  a  sample  step  function  g(t),  and  its  linear  approxi- 
mation h(t)a   Using  h(t),  the  generic  integral,  N,  becomes 
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TCON(l)   _       /RTCON(l)   ....        /ALTIM 

s°dt  ♦   /         sMt)dt  ♦   /      sh(t)dt 
JTCON(I)  JRTCON(I) 

and,  by  equation  (307),  becomes 

H(4)   H(3)     TRrniv.    H(6)   H(5) 

.*,.*«,.  /,%     -rr-^r,.,  s       -s  TBCON     s       -  3 

=  TC0N(1)  +  TBCON — r-r +  — —   7— T"T 

Ln ( s )  2       Ln  ( s ) 

The  following  derivatives  are  of  immediate  interest  in  evaluating  the  inte- 
grals Pr(A), 

/   m  \             Tprnw   i  *  l  ruf^i   H(4)-l   u/,,    H(3)-l.    ,  H(4)-l  H(3)-l. 
f   N   1     TBCON   Ln ( s )  ( H ( 4 )  s -  HO)   s )  -  (s ^s ) 

ALTIM   =  ALTIM  .  2,  > 

\  I  Ln  (s) 

Tormi   .  I    \    (u(c\      H(6)-l   „,-,   H(5)-l.    ,  H(6)-l  H(5)-l. 
TBCON   Ln(s)  (H(6)  s -  H(5)  s )  -  (s -s ) 

+  2  ALTIM  ,  2,  . 

Ln  (s) 


and 


D 


)ALTIM-N     \  ^  ALTIM-  U°S)  D  IaZtIhI 

ALTIM  (1-s)  ,,    .2 

'  \i   -  s; 


In  terms  of  these  values,  the  integrals  Pr(A)  become 

(3.9)      PA1  =  Pr  (a  submarine  initially  on  station  survives  the  first 

Markov  chain  to  attempt  passage  home) 


N  ALTIM  -  N 


*  P 


ALTIM     K5    1-s 

.        .  Q 

(3.10)     PA2  =  Pr  (a  submarine  initially  on  station  fires  one  salvo  at 

the  cargo  ships  in  the  first  Markov  chain) 

.       .    ALTIM  -  N 
=   lp3  +  V    ALTIM  (1  -  s) 
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TCOM(i) 


RTCON(l) 


ALTIM 


FIGURE  4 
A  SAMPLE  STEP  FUNCTION,  g(t),  AND  ITS  APPROXIMATION,  h(t) 
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d 
(3.11)     PA3  =  Pr  (a  submarine  initially  on  station  fires  two  salvos  at 

the  cargo  ships  in  the  first  Markov  chain) 


,       ,     ALTIM  -  N 
=   (PK  +  Pj 


5   p6'   ALTIM  (1  -  s) 

(3.12)  PA4  =  Pr  (a  submarine  initially  on  station  survives  the  inter- 

action) 

t-,  n\     N      ,  n.     ALTIM  -  N 

=   U  +  P3  D)   ALTIM  *  (p5  +  P3  q3  D)    ALTIM(1  -  s) 

(3.13)  PA5  =  Pr  (a  submarine  initially  on  station  fires  on  cargo  ships 

in  the  second  Markov  chain) 

(  ^    n     ALTIM  -  N 

=   (q3  *  V  P3  °   ALTIM  (1  -  s) 

This  completes  the  evaluation  of  the  integrals  Pr(A)  for  those  submarines 
initially  on  station. 

The  integrals  Pr(A)  are  trivially  Stieltjes  integrals  for  those  sub- 
marines which  arrive  on  station  after  the  war  starts.   The  following  is  based 
on  the  assumption  that  the  arrival  time  for  these  submarines  is  later  than 
RTCON(l),  i.e.,  the  full  convoy  arrival  rate  of  two  convoys  per  TBCON  hours 
has  been  achieved  prior  to  submarine  arrival  at  the  subline.   In  this  case, 

Pr(A)   =  Pr(A/  g (maximum  value  of  t)) 
but 

g(max  t)   = 


ALTIM  -  submarine  transit  time  to  station 

x  2 


TBCON 


Neglecting  possible  delays  engendered  during  HUK  attacks  and  during  the  barrier 

D5UBSTA 
penetration,  the  submarine  transit  time  is   ciracD. —  ,  therefore 
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Pr(A)   =  Pr(A/  n+1  = 


ALTIM 


D5UB5TA 
SUBSPD 


TBCON 


x  2 


Letting  KONE  = 


ALTIM  - 


DSUBSTA 
SUBSPD 


TBCON 


x  2 


,  the  following  probabilities  can 


be  found 

(3.14)     PA6  =  Pr  (a  submarine  not  initially  on  station  survives  the  first 

Markov  chain  to  attempt  passage  home) 

KONE       1  -  sK0NE 
=  s     +  pc 


5  1  -  s 


(3.15) 


(3.16) 


(3.17) 


pa7  i 


Pr  (a  submarine  not  initially  on  station  fires  one  salvo 
at  the  cargo  ships  in  the  first  Markov  chain) 

KONE 


(P3  ♦  P4) 


1-3 
1  -  S 


PA8  = 


Pr  (a  submarine  not  initially  on  station  fires  two  salvos 
at  the  cargo  ships  in  the  first  Markov  chain) 
KONE 


(P5  +  P6) 


1-3 
1-3 


d 
PA9  = 


Pr  (a  submarine  not  initially  on  station  survives  the  sub- 
marine-convoy interaction^) 

,nuc  .    KONE 

=   (1  ♦  p3  D)   s      ♦  (p5  +  p3  q3  D)   — - 

(3.18)    PA10  =  Pr  (a  submarine  not  initially  on  station  fires  on  cargo 

ships  in  the  second  Markov  chain) 
KONE 

■  (q3  *  V  p3  D  r^f 


This  completes  the  evaluation  of  the  integrals  in  terms  of  the  evaluated  or 
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known  FORTRAN  variables  and  the  as  yet  unevaluated  transition  probabilities. 

3o2„2  The  Barrier  Penetration,  Transition  Probabilities  and  Related  Expected 
Numbers 

Until  now,  it  has  been  convenient  to  suppress  in  the  analysis  any  nota- 
tion concerning  submarine  type  and  sonar  condition.   The  Markov  chain  notation 
was  difficult  enough  without  the  burden  of  two  additional  subscripts.   The 
remainder  of  this  analysis  will  distinguish  submarines  by  type  and  by  assigned 
sonar  condition.   It  will  be  seen  that  each  of  the  probabilities  given  by 
equations  (3.9)  through  (3.18)  is  a  function  of  submarine  type  and  sonar 
condition. 

Placement  of  submarines.   In  SEALIFT,  there  are  eight  ways  to  classify 
submarines  of  a  particular  type  by  placement  at  the  start  of  war:   two  choices 
for  on  station  or  not  and  four  choices  of  sonar  box.   This  classification  is 
done  by  Monte-Carlo  methods.   Type  sonar  condition  within  each  box  is  de- 
termined by  input.   For  the  purposes  of  this  analysis,  sonar  box  classifica- 
tion is  immaterial  except  for  the  restriction  that  the  fourth  sonar  box  is 
not  included  in  the  analysis.   For  the  initial  placement  of  submarines  the 
expected  number  of  type  I  submarines  initially  on  station  in  type  K  sonar 
conditions  is  given  by 

E2(I,K)   =  NSUB'I'T  x  FSUB'I'OS  x  ^P'l'BOXU) 
where  the  summation  is  over  all  sonar  boxes  assigned  sonar  condition  K.   The 
expected  number  of  type  I  submarines  not  initially  on  station  but  assigned 
to  type  K  sonar  conditions,  E3(I,K),  is  derived  from  two  sources,  the  fraction 
of  on  hand  submarines  not  initially  on  station  and  the  submarines  built  during 
the  war.   Using  a  nominal  thirty  day  month, 

E3(I,K)   =   (NSUB'I'T  (1  -  FSUB'I'OS)  +NSUB'!'  x  ^j^-)  x  ]Tp«I'BOX(J) 
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where  the  summation  is  over  all  sonar  boxes  assigned  sonar  condition  K» 

Barrier  penetration,,   A  submarine  will  successfully  transit  the  eight 
distinct  barriers  if  and  only  if  she  successfully  transits  each  barrier.   In 
the  5EALIFT  game,  the  barriers  are  evenly  divided  between  two  types,  constant 
and  variable.   The  latter  type  are  variable  in  the  sense  that  'the  barrier- 
detects-submarine  probability  will  vary  with  barrier  strength.   Under  the 
basic  assumptions  used  in  this  expected  value  analysis,  however,  each  of 
the  barriers  is  treated  as  a  constant  barrier.   For  constant  barriers  there- 
fore 

PA11(I,J)   =  Pr  (a  type  I  submarine  successfully  transits  the  J  th 
constant  barrier) 
=  1  -  Pr  (she  is  sunk  in  this  barrier) 
=  1  -  PBARSC(I,J) 

and 

.      d 
PA12(I,J)   =  Pr  (a  type  I  submarine  successfully  transits  the  J  th 

variable  barrier) 

=  1  -  Pr  (she  is  sink  in  this  barrier) 

=  1  -  PBARSV(I,J)  x  (PVBDET(I,J)  +  VAIRDET(J)) 

Therefore  the  probability  that  a  type  I  submarine  successfully  transits  the 

barriers,  P0(I),  is 

4 

P0(I)    =  Tf      PA11(I,J)  x  PA12(I,J) 
J=l 

Convoy  detection.   In  the  SEALIFT  game,  all  convoys  cross  the  subline 
at  the  subline' s  midpoint a   Detection  probabilities  are  of  the  "cookie 
cutter"  type  where  submarines  are  placed  (at  each  convoy  crossing)  uniformly 
on  the  subline.   The  probability  an  (I,K)  submarine  detects  a  convoy  on  any 
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one  trial,  1  -  P1(I,K),  is  then 

,    ,     2  x  RGDT(K.I) 
1  -  piu,k;  -    5UBLINE 

The  transition  probabilities,  p   =  q   =  P1(I,K),  are  thereby  fully  defined, 
and  are  dependent  on  submarine  type  and  on  sonar  condition. 

Penetration  of  the  forward  screen  escort.   This  event  and  the  following 
events  are  the  substance  of  Event  39  (Submarine  versus  convoy)  of  the  SEALIFT 
game.   Figure  5  is  the  simplified  truth  table  (tree)  of  that  part  of  Event  19 
from  which  the  F5E  associated  probabilities  and  expectations  used  in  this 
analysis  are  derived,,   All  of  these  probabilities  are  conditioned  on  the 
submarine  having  detected  a  convoy.   The  input  parameters  FAIR5C  and  ATSFE(I) 
are  obviously  not  probabilities.   FAIRSC  is  transformed  into  a  0,1  parameter 
in  this  analysis  by  the  transformation  2  -  FAIRSC,  allowing  its  use  in  prob- 
ability statements.   ATSFE(I)  is  another  tactical  limitation  of  the  type  I 
submarine  and  this  parameter  is  usually  a  "small"  integer  for  the  conventional 
submarine  and  a  "large"  integer  for  the  nuclear  submarine.   If  ATSFE(I)  =  0, 
a  type  I  submarine  will  not  be  allowed  to  attempt  repenetration  of  the  FSE 
in  the  event  she  is  rebuffed.   For  this  submarine  in  type  K  sonar  conditions 
the  following  probabilities  are  obtained 

P7(I,K)   =  Pr  (a  submarine  penetrates  the  FSE;  ATSFE(I)  =  0) 

=  1  -  (PDFS(I,K)  +  PDFER) 
P9(I,K)   =  Pr  (a  submarine  returns  to  station;  AT5FE ( I )  =  0) 

=   (PDFS(I,K)  +  PDFER)  x  (1  -  PDFCE(I,K))  x  (1  -  (2  -  FAIRSC) 
x  PSFSE(I)  x  PFCES(I.K)  x  PSFCE(I,K)) 


and 


P8(I,K)   =  Pr  (a  submarine  is  killed  attempting  FSE  penetration;  AT5FE(I)»0) 
=  1  -  P9(I,K)  -  P7(I,K) 
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IS  SUBMARINE  DETECTED 
BY  FSE? 


YES 

PDFS(I,K) 
+  PDFER 


NO 


© 


IS  SUB  SUNK? 


YES  NO 

PDFCE(I,K) 


© 


DOES  SUB  ATTACK  FSE? 


YES 

(2  -  FAIRSC) 

x  PSFSE(I) 


NO 


DOES  SUB  SINK 
AN  FSE  SHIP? 


YES        NO 
PFCES(I,K) 


IS  SUB  SUNK 
AFTER  FLAMING 
DATUM? 


YES       NO 
PSFCE(IPK) 


© 


Legend: 
A:  submarine 
B:  submarine 
X:   submarine 


CAN  SUB  ATTEMPT  REPENETRATION? 


NO 


YES 


NUMBER  OF  ATTEMPTS  IS 
A  PRE-SET  DECISION  GIVEN 
BY  ATSFE(I) 


penetrates  FSE 
returns  to  station 
is  killed 


FIGURE  5 
ARBORESCENCE  ILLUSTRATING  THE  SUBMARINE-FORWARD  SCREEN  ESCORT  INTERACTION 
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In  general,  if  ATSFE ( I )  =  n,  n  attempts  at  penetration  will  be  allowed,  and 
P7(I,K,n)   =  Pr  (a  submarine  penetrates  the  F5E ;  ATSFE(I)  =  n) 


_  P7(I  K)   1  -  (P9(ItK)) 
"  P7U'K)    i  -  P9(I,K) 


n+1 


and 


P9(I,K,n)   =  Pr  (a  submarine  returns  to  station;  ATSFE(I)  =  n) 
=   (P9(I,K))n+1 


P8(I,K,n)   =   Pr  (a  submarine  is  killed  attempting  F5E  penetration; 


ATSFE(I)  =  n) 

n+1 


pnn  K)       i  -  (P9U.K)) 

-   PB(I.K)     1  _  p9(I)(<) 


As  derived  herein  these  values  are  predicated  only  on  the  approximation  that 
the  F5E-detects-submarine  probability  has  not  been  degraded  by  prior  sub- 
marine attack  either  by  the  submarine  under  consideration  (ATSFE(I)  >  0)  or 
by  earlier  submarines  and  that  when  the  submarine  attempts  to  attack  the 
screen  she  will  have  sufficient  torpedoes  to  do  so.   To  find  the  expected 
number  of  torpedoes  a  submarine  fires  at  FSE  ships  given  she  penetrates  the 
screen,  suppose  ATSFE (I)  =  n,  and  define  the  following  events 

C  .   =  event  that  a  submarine  penetrates  the  screen  on  the  i  th 
J 

try,  j  =  1,  .oo ,  n+1 
D.   =  event  that  a  submarine  fires  i  salvos,  each  of  size  ASALV0(I! 

during  an  attempted  penetration,  i  =  0,  ...,  n+1 
Then  E1(I,K)  can  be  defined  as  the  expected  number  of  torpedoes  fired  by  an 
(I,K)  submarine  during  an  attempted  penetration  of  the  FSE  given  she  penetrates 
the  FSE  and  ATSFE (I)  =  n,  i0e0, 
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(3.19)  n+1 

n        L         Pr(C,  D.) 
E1(I,K)   =   ASALVO(I)    )    i   ifili. 


i=0         P7(I,K,n) 

The  probability  Pr  (C.,  D.)  can  be  derived  by  the  following  argument:   the 

event  (C.,  D.)  is  a  sequence  of  j  trials  where  the  j  th  and  last  trial  is 
J   i 

a  non-detection  by  the  FSE.   This  trial  is  preceeded  by  i  trials  (not 
necessarily  in  sequence)  in  which  the  submarine  is  detected,  not  sunk,  does 
attack  the  FSE ,  and  either  sinks  an  F5E  ship  and  is  not  killed  or  does  not 
sink  an  FSE  ship.   The  last  trial  is  also  preceeded  by  j  -  i  -  1  trials  in 
which  the  submarine  is  detected,  not  sunk,  and  does  not  attack  the  FSE. 
Denote  the  probability  of  each  type  trial,  in  the  order  listed,  by  Tl,  T2, 
and  T3.   Then 

T1(I,K)   =  P7(I,K) 

T2(I,K)   =   (PDFS(I,K)  +  PDFER)  x  (1  -  PDFCE(I,K))  x  (2  -  FAIRSC) 
x  P5FSE(I)  x  (1  -  PFCES(I,K)  x  PSFCE(I.K)) 
and 

T3(I,K)   =   (PDFS(I,K)  +  PDFER)  x  (1  -  PDFCE(I,K))  x  (1  -  (2  -  FAIRSC) 
x  PSFSE(I)) 
Therefore, 

Pr(C,  D.)   =   (J_1)  x  (T2(I,K))i  x  (T3(  I  ,K) )  j"i"1  x  P7(I,K) 

and   equation    (3.19)    becomes 

,  n+1 

(3°19)  n  I         M:1)    (T2(I,K))ix    (TSd.K))^1"1 


E1(I,K)      =     ASALVO(I) 


I   i 


■i=i+l 


1 


i=l  1  -   (P9(I,K))rH'1 

1  -  P9(I,K) 
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Penetration  of  the  inner  screen  (ASW).   Having  penetrated  the  FSE ,  an 
attacking  submarine  must  penetrate  the  inner  screen  prior  to  reaching  the 
cargo  ships „   Figure  6  is  the  simplified  tree  of  that  part  of  Event  19 
(Submarine  versus  convoy)  of  5EALIFT  from  which  the  ASW  associated  probabil- 
ity of  penetration  is  derived  for  this  expected  value  analysis.   This  portion 
of  the  analysis  is  based  upon  an  undated  and  informal  modification  of  Event 
19  which,  however,  modifies  only  the  ASW  penetration  portion.   The  assump- 
tions upon  which  the  tree  is  constructed  are  the  same  as  used  for  the  FSE 
penetration  analysis.   That  is,  the  ASW-detects-submarine  probability  is 
assumed  not  degraded  by  prior  ASW  ship  losses  and  that  if  a  submarine  attacks 
the  ASW,  she  is  assumed  to  have  sufficient  torpedoes  to  do  so.   Because  of 
the  method  used  in  deriving  the  ASW  penetration  probability  for  this  analysis, 
parasitic  torpedo  usage  could  not  be  determined.   It  was  observed  that  two 
mutually  exclusive  events  prevailed,  either  the  submarine  penetrated  the  ASW 
or  she  was  sunk.   The  method  used  for  computing  the  probability  of  these 
events  occurring  was  to  first  determine  the  probability  the  submarine  was 
sunk.   The  desired  probability  of  penetration  is  found  by  complementation. 

The  various  paths  to  "submarine  death"  during  the  ASW  penetration  are 

stylized  in  Figure  7.   The  circled  integers,  (l)»  and  the  circled  arabic  A, 

(a),  in  Figure  7  represent  branch  points  corresponding  to  those  of  the  tree, 

Figure  6.   The  circled  arabic  X,  (x),  represents  the  event  "submarine  is 

killed."  The  symbol  p   is  the  probability  of  transition  from  branch  point 

r  to  branch  point  s,  indexed  by  the  integer  i.   From  Figures  6  and  7  the 

probabilities  for  an  (I,K)  submarine  are  read  off; 

Plil   =  PDETU(I,K)  x  VA(K) 
r14 

PA8  =  PBB(I,K)  x  PAT(I,K) 
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IS  SUB  DETECTED  BY  BOTTOM  BOUNCE? 


YES 
PBBdTk) 


NO 


DOES  SUB  ATTACK  UNALERTED  INNER 
SCREEN? 


YES 
PSUBAT(I) 


NO 


IS  A  SCREEN 
SHIP  SUNK? 


YES 
PKOSC(I) 


NO 


DOES  HIGHLY  ALERTED 
INNER  SCREEN  DETECT 
SUB? 


YES 

PDETU(I,K) 
x  VM(K) 


NO 


(  MR«*-1  ) 


© 


DOES  UNALERTED 
INNER  SCREEN 
DETECT  SUB? 


YES        NO 
PDETU(I,K) 

(  MR^l  ) 


0 


QTU(389)  =  1? 


YES 


NO 


(  MR-»-l  ) 


DOES  ASW  SCREEN 
ATTACK  SUB? 


YES 
PAT(I,K) 


NO 


© 


© 


DOES  ALERTED  SCREEN 
DETECT  SUBMARINE? 


YES 

PDETU(I,K) 
x  VA(K) 


NO 


(  MR--1  ) 


Legends 

J:   Submarine  penetrates 

ASW 
Xs   Submarine  is  killed 

INTEGER:   Branch  point 


Q 


FIGURE  6 
ARBORESCENCE  ILLUSTRATING  THE  SUBMARINE-ASW  SCREEN  INTERACTION 
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8  >*■ 


IS  SUB  KILLED? 


YES 

PKBB(I,K) 


NO 


© 


DOES  SUB  ATTACK  THIS 
GROUP? 


YES 
PABB(I,K) 


NO 


DOES  SUB  KILL  A 
SCREEN  SHIP? 


YES        NO 
PKBBG(I) 


DOES  THE  SCREEN  COUNTER  ATTACK? 


YES 
PCC(I,BB) 


NO 


BB  — BB  +  1 


BB  >  3? 


YES    NO 


© 


FIGURE  6  (continued) 
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DOES  SCREEN  KILL  SUB? 

YES                  NO 
PKINS(I) 

© 

' 

DOES  SUB  ATTACK  SCREEN? 

YES              NO 
PAA(I) 

' 

DOES  SUB  SINK  A 
SCREEN  SHIP? 

YES       NO 
PKSC(I) 

_J 

QTCU77)  =  1? 

NO      YES 

* 

1 

\ 

y^j 

DOES  SCREEN  COUNTER  ATTACK? 

YES                     NO 
PSCA(I,MR) 

\ 

Q 

D 

(  MR --MR  +  1  ) 


MR  >  3? 


YES    NO 


© 


FIGURE  6  (continued) 
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© 


DOES  ALERTED  SCREEN 
DETECT  SUB? 


YES 

PDETU(I,K) 
x  VA(K) 


NO 


(  MR— -1  ) 


© 


FIGURE  6  (continued) 
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FIGURE  7 
STYLIZED  PATHS  TO  SUBMARINE  DEATH  DURING  THE  SUBMARINE-ASW  SCREEN  INTERACTION 
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p8x  =  PKBB(I,K) 

p*   =   (1  -  PKBB(I,K)J  x  PCC(I.i)  x  (1  -  PABB(I,K)  +  PABB(I,K)  x 


(PKBBG(I')  +  (1  -  PKBBG(I))  x  (QTC(177) ) ) ) ,   i=l 3 


PB1  =  X  -   PB8  "  P8X'   i=1'  •••»  3 

P81  =  l  "   PKBB(I'K) 

p*  =  (1  -  PKINS(I))  x  PSCA(I,i)  x  (1  -  PAA(I)  x  (PK5C(I) 
r44 


+  (1  -  PK5C(I))  x  (1  -  QTC(177)))),   i=l,  ...,  3 


p4x  =  PKINS(I) 


P    =   (1  -  PBB(I,K))  x  PSUBAT(I)  x  (1  -  PKOSC(I))  x  (1  -  QTU(389)) 

MX 


and 


P    =  PBB(I,K)  x  (1  -  PAT(I,K))  x  PDETU(I,K)  x  VA(K) 
A4 

+  (1  -  PBB(I,K))  x  (PSUBAT(I)  x  (PKOSC(I)  x  PDETU(I,K)  x  VM(K) 

+  (1  -  PKOSC(I))  x  QTU(389))  +  (1  -  PSUBAT(I))  x  PDETU(I,K))) 

The  probability  that  the  submarine  penetrates  the  inner  screen,  P3(I,K),  can 
then  be  found  as 

P3(I,K)   =  1  -  Pr  (Submarine  is  sunk  attempting  A5W  penetration) 
=  1  -  Pr  (a  transit  from  the  first  f4j  to  (x)  by  all 

possible  paths)  x  Pr  (a  transit  from  TaJ  to  first 
\Aj   by  all  possible  paths)  +  Pr  (a  transit  from 
(a)  to  fXJ  by  all  possible  paths  not  passing 
through  any  [A J   ) 
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Each  of  these  probabilities  can  be  read  from  Figure  8  to  be  the  following: 
Pr  (a  transit  from  the  first  f4j  to  (XJ   by  all  possible  paths) 

=  P4X  +  P44  (P4X  +  P44  (p4X  +  P44  P4XJ) 
Pr  (a  transit  from  (a)  to  first  (Aj  by  all  possible  paths) 

=  PA1  P14  +  PAB(P81  P14  +  P8B(p81  P14  *  P88lp81  P14  +  P88  P81  P14})) 

and 

Pr  (a  transit  from  (a j  to  fxj  by  all  possible  paths  not  passing  through 
any  (4)  ) 

=  PA8  (P8X  +  P88(P8X  +  P88(p8X  +  P88  P8X)}) 

Convoy  attack  phase.   It  is  particularly  in  the  convoy  attack  that  the 
flexibility  of  the  SEALIFT  game  complicates  the  expected  value  analysis. 
Possible  realizations  of  this  phase  depend  heavily  upon  the  chosen  war 
option  as  listed  in  SEALIFT  I  1  and  upon  the  various  torpedo  salvo  sizes 
chosen  for  each  submarine  type,,   In  view  of  these  considerations  and  the 
problem  which  prompted  this  expected  value  analysis,  a  particular  character- 
ization of  the  convoy  attack  phase  has  been  chosen.   It  is  assumed  that 

1.  one  type  of  torpedo  will  be  used  against  both  cargo  and 
screen  ships, 

2.  the  salvo  size  against  cargo  ships  is  C5ALV0(I)  and  that 
CSALVO(I)  =  -^  the  initial  torpedo  supply, 

3„   the  salvo  size  against  screen  ships  is  ASALVO(I)  and  that 

ASALVO(I)  1  1/6  x  CSALVO(I), 
4.   in  firing   one  CSALVO  against  cargo  ships,  the  combination 

of  CSALVO(I)  and  the  probability  of  kill  per  torpedo 
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essentially  eliminates  the  probability  of  sinking  no  cargo 
ships,  and 
5.   having  fired  two  salvos  at  cargo  ships,  the  submarine's 
torpedo  supply  is  exhausted. 
Based  on  these  assumptions,  Figure  8  is  the  simplified  tree  of  that 
part  of  Event  19  (Submarine  versus  convoy)  of  SEALIFT  from  which  the  convoy 
attack  associated  probabilities  are  derived  from  the  first  Markov  chain  in 
this  analysis.   For  an  (I,K)  submarine  sonar  conditions  the  conditional 
probabilities  of  interest  are 

P41(I,K)   =  Pr  (a  submarine  shoots  one  salvo  at  cargo  ships  in  the 
first  Markov  chain) 
=  PD5AU(I,K) 
P42(I,K)   =  Pr  (a  submarine  shoots  two  salvos  at  cargo  ships  in  the 
first  Markov  chain) 
=  1  -  PDSAU(I,K) 
P5(I,K)   =  Pr  (a  submarine  escapes/  she  shoots  one  salvo  at  cargo 
ships  in  the  first  Markov  chain) 
=  1  -  PSSAU(I,K) 
and 

P6(I,K)   =  Pr  (a  submarine  escapes/  she  shoots  two  salvos  at  cargo 
ships  in  the  first  Markov  chain) 
=  1  -  PDSAU(I,K)  x  PSSAU(I,K) 
If  the  expected  number  of  torpedoes  fired  by  an  (1,K)  submarine  at 
cargo  ships  in  the  first  Markov  chain  given  the  submarine  fires  one  and  two 
salvos  are  E4(I,K)  and  E5(I,K)  respectively,  then  E4  and  E5  are  roughly 
approximated  by: 

E4(I,K)   t  CSALVO(I) 

49 


SUB  SHOOTS  AT  CARGO  SHIPS 


IS  SUB  DETECTED  BY  SAU? 


YES 
PD5AU(I,K) 


IS  SUB  SUNK? 


YES         NO 
P5SAU(I,K) 


0 


DOES  SUBMARINE 
COUNTER  ATTACK? 


YES         NO 

DECISION 
ATTSAU(I)  =  1 


SUB  SINK  A  SAU 
SHIP? 


YES         NO 
PSAUS(I,K) 


SUB  REATTACK  CONVOY? 


NO  YES 

NUMBER  OF  ATTEMPTS  IS 
A  PRE -SET  DECISION 
GIVEN  BY  RA5W(I) 


SUB  ESCAPES 


NO 


SUB  SHOOTS  AT  CARGO  SHIPS 


IS  SUB  DETECTED  BY  SAU? 


NO 


YES 
PDSAU(I,K) 


SUB  ESCAPES 


SUB  SUNK? 


YES      NO 
PSSAUS(I,K) 


© 


SUB  ESCAPES 


NOT  CONSIDERED 


FIGURE  8 
ARBORESCENCE  ILLUSTRATING  THE  CARGO  SHIP  ATTACK  PHASE 
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and 


E5(I,K)   t  TORPS(I,ASW)  -  CSALVO(I)  -  E1(I,K) 


The  tree  associated  with  the  cargo  attack  phase  for  those  submarines 
which  successfully  enter  the  second  Markov  chain  is  obtained  from  Figure  9 
by  changing  the  entire  right  branch  at  the  branch  point  K  into  the  event 
"Submarine  Escapes."  The  conditional  probabilities  of  interest  are  then 

P4(I,K)   =  Pr  (a  submarine  shoots  her  remaining  torpedoes  at  cargo 
ships) 
=  1 
Pll(IfK)   ss  Pr  (a  submarine  escapes  given  she  shoots) 
=  P5(I,K) 
The  expected  number  of  torpedoes  a  submarine  fires  at  cargo  ships  in 
the  second  Markov  chain  given  she  fires,  E6(I,K),  is  roughly  approximated  by 
E6(I,K)   I  T0RP5(I,ASW)  -  CSALVO(I)  -  2  x  E1(I,K) 
The  transition  probabilities.   The  calculation  of  the  Markov  chain 
transition  probabilities  needed  for  this  expected  value  analysis  can  now  be 
obtained  by  combining  the  conditional  probabilities  derived  in  the  previous 
sections.   These  transition  probabilities  are: 


Pl  =  qi  =  P1(I,K) 


p2   =  q2   =   (1  -  P1(I,K))  x  P9(I,K,n) 


p3  =   (1  -  P1(I,K))  x  P7(I,K,n)  x  P3(I,K)  x  P41(I,K)  x  P5(I,K) 


p4  =   (1  -  P1(I,K))  x  P7(I,Ktn)  x  P3(I,K)  x  (1  -  P5(I,K)) 

pc   =   (1  -  P1(I,K))  x  P7(I,K,n)  x  P3(IfK)  x  P42(I,K)  x  P6(I,K) 

p6  =   (1  -  P1(I,K))  x  P7(I,K,n)  x  P3(I,K)  x  P42(I,K)  x  (1  -  P6(I,K)) 
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p?  =  q5  =   (1  -  P1(I,K))  x  P7(I,K,n)  x  (1  -  P3(I,K)) 

pQ  =  q6   =   (1  -  P1(I,K))  x  P8(I,K,n) 

q3  =  (1  -  P1(I,K))  x  P7(I,K,n)  x  P3(I,K)  x  P4(I,K)  x  P5(I,K) 

q4  =  (1  -  P1(I,K))  x  P7(I,K,n)  x  P3(I,K)  x  P4(I,K)  x  (1  -  P5(I,K)) 

3.2.3  Expected  Number  of  Torpedoes  Fired  at  Cargo  Ships 

In  order  to  determine  the  expected  number  of  torpedoes  fired  at  cargo 
ships  for  any  given  cycle,  define  the  following  to  be 

F(u,v,I,K)   =  event  an  (I,K)  submarine  fires  u  salvos  at  cargo  ships 
in  the  v  th  Markov  chain,   u,  v  =  1,  2 
NS(w,I,K)   =  the  number  of  (I,K)  submarines  which  arrive  on  station 

for  the  w  th  cycle,   w  =  1,  2,  3,  . .. 
NT(w,I,K)   =  the  total  number  of  torpedoes  fired  at  cargo  ships  during 
the  w  th  cycle  by  (I,K)  submarines,   w  =  2,  3,  ... 
The  following  equation  is  then  used 
(3.20) 

E(NT(w,I,K))   =  E(NS(w,I,K))  x  ( Pr(F(l,l,I ,K) )  x  E4(I,K) 

+  Pr(F(2,l,I,K))  x  E5(I,K) 
+  Pr(F(l,2,I,K))  x  E6(I,K)j 
For  the  first  cycle,  distinction  is  made  for  the  initially  on-station, 
initially  not-on-station  contingency.   For  each  submarine  type  and  for  each 
sonar  condition  the  expected  number  of  torpedoes  fired  in  each  contingency 
is  computed  as 

E (NT. (1, I ,K) )   =  E  (number  of  torpedoes  fired  by  (I,K)  submarine  initially 

on  station  during  the  first  cycle) 
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=  E2(I,K)  x  (  PA2(I,K)  x  E4(I,K) 
+  PA3(I,K)  x  E5(I,K) 
+  PA5(I,K)  x  E6(I,K)J 
and 

E(NT_(1,I ,K) )   =  E  (number  of  torpedoes  fired  by  (I,K)  submarines 

initially  not-on-station  during  the  first  cycle) 
=  E3(I,K)  x  P0(I,K)  x  j  PA7(I,K)  x  E4(I,K) 

+  PA8(I,K)  x  E5(I,K) 
+  PA10(I,K)  x  E6(I,K)( 

The  expected  number  of  survivors  from  the  first  cycle,  i«e.,  those  who  safely 
reach  home  port,  EC(1,I,K),  is 

EC(1,I,K)   =   (E2(I,K)  x  PA4(I,K)  +  E3(I,K)  x  PO(I,K)  x  PA9(I,K))x  P0(I,K) 

According  to  the  rules  of  SEALIFT,  the  survivors  are  re-apportioned  among 
the  sonar  boxes  to  maintain  balance „   The  expected  number  of  re-apportioned 
survivors  of  the  (I,K)  classification,  EC(2,I,K),  is 


EC(2,I,K)   =  EC(1,I,K)  x  Vp'I'BOX(J) 


where  the  sum  is  over  those  sonar  boxes  of  like  sonar  conditions „   To  deter- 
mine E(NT(2,I  ,K) ) ,  equation  (3.20)  is  applied  again,  and  becomes 

E(NT(2,I,K))   =  ^yj1^  x  E(NT2(1,I,K)) 
The  expected  number  of  survivors  from  the  second  cycle,  EC(3,I,K),  is 
EC(3,I,K)   =  EC(2,I,K)  x  PA9(I,K)  x  P0(I,K)  x  P0(I,K) 
For  the  w  th  cycle,  w  =3,  the  expected  number  of  re-apportioned  survivors 
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from  the  previous  cycle,    EC(2w-2, I ,K) ,    serve   as  the   submarine   input,    and 
the   following   iteration   formulae   are   found, 


EC(2w-2,I,K)      -     EC(2w-3,I,K)    »  §{§%$ 


and 


E(NT(w,I,K))      = 


rr,,      o    t    v\         E(NT(w-l,I,K)) 
EC(2w-2,I,K)    x     EC(2w,4>I>K) 


The  expected  number  of  survivors  from  the  w  th  cycle  (not  yet  re-apportioned 
among  the  sonar  boxes),  EC(2w-l,I,K)  is 


EC(2w-l,I,K)   =  £C(2w-2,I,K)  x 


EC(2w-3.I.K) 
EC(2w-4,I,K) 


An  estimate  of  the  number  of  cycles  for  the  war  is  given  by  the  minimum 
number  of  cycles  a  live  submarine,  not  initially  on  station,  might  expect  to 
achieve o   This  number  is  a  function  of  submarine  type  and  is  given  by  the 
integer  L(I) 


L(I) 


TWAR 


ALTIM(I)  +  DTSPRT(I)  + 


D5UB5TA 
SUBSPD(I)  _ 


An  estimate  for  total  expected  number  of  torpedoes  fired  at  cargo 
ships  by  (I,K)  submarines  for  the  war,  E (NT (I  ,K) ) ,  is  then 

f 

E(NT(I,K))   =  E(NT  (1,I,K))  +  E(NT,(1, I ,K) )  +  /    E(NT(w,I,K) 

1  2  *—*> 

w=2 

Multiplying  this  number  by  the  single  shot  kill  probability,  PSCS(I,K),  pro- 
vides an  estimate  for  the  number  of  cargo  ships  sunk  by  (I,K)  submarines, 
E(CS(I,K)), 

E(CS(I,K))   =  E(NT(I,K))  x  PSCS(I,K) 
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Finally,  summing  E(CS(I,K))  over  the  indicies  I  and  K  provides  an  estimate 

for  the  total  expected  number  of  cargo  ships  sunk  during  the  entire  war,  E(CS) 

3    2 

E(CS)   =   }_   2_   E(CS(I,K)). 
1*1  K=l 
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4.   Conclusions  and  Recommendations 

As  a  model  purporting  to  represent  a  probabilistic  real  world  situation, 
the  5EALIFT  game,  the  analysis  presented  in  this  thesis  remains  largely  un- 
tested.,  The  opportunity  and  data  base  for  testing  the  model's  reliability, 
accuracy  and  sensitivity  to  the  proposed  approximations  are,  however,  readily 
available,,   If  the  5EALIFT  game  is  contemplated  as  an  analytical  tool  in  the 
CYCLOPS  III  study,  the  accuracy  of  this  expected  value  model  should  be  ex- 
plored.  An  acceptably  reliable  programmed  version  of  this  model  should 
effect  appreciable  savings  in  computer  timec 

For  her  continued  patience,  encouragement  and  good  humor,  I  acknowledge 
the  support  of  my  wife  in  the  preparation  of  this  thesis. 
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APPENDIX  I 

DETAILED  DESCRIPTION  OF  THE  5EALIFT  GAME 

The  following  description  is  quoted  from  SEALIFT  I  ll   : 

II.   GENERAL  DESCRIPTION 

Convoys 

In  this  game,  friendly  forces  have  two  ports  to  transit  between  home 
port  and  delivery  porto   During  this  transit  the  convoy  may  cross  up  to  two 
"constant'*  barriers  (oneway  transit).   The  constant  barrier  in  this  simulation 
can  sink  ships  that  cross  it  but  cannot  itself  be  weakened  by  these  sinkings 
or  by  the  actions  of  the  ships  crossing  it.   Mine  fields,  air  barriers,  etc., 
can  be  described  as  "constant"  barriers.   The  "variable"  barrier,  used  against 
submarines,  can  sink  ships  that  traverse  it  but  can  be  weakened  by  having  the 
ships  that  comprise  it  sunk.   Barriers  of  this  type  might  be  submarine  bar- 
riers, surface  ship  barriers,  etc.   The  "variable"  type  of  barrier  will  be 
described  in  greater  detail  later  in  this  section.   It  should  be  pointed  out 
that  while  each  barrier  that  could  be  encountered  by  a  convoy  or  submarine 
will  be  described  herein,  any  one  or  all  of  these  barriers  may  be  removed  for 
a  desired  play  simply  by  setting  to  zero  its  probability  of  detecting  a  pro- 
spective "transitor". . . 

At  certain  times,  which  are  inputs,  a  convoy  is  formed  in  the  home  port. 
The  convoy  is  made  up  of  three  types  of  ships:   cargo  ships,  ASW  escort  ships, 
and  forward  screen  ships*.   Each  convoy  when  it  forms  draws  its  required 
complement  of  ships  from  a  "port  pool."  This  pool  is  stocked  by  ships  from 
returning  convoys  and  by  the  ship  building  rate.   If  the  required  ships  are 
available,  the  convoy  sails  toward  its  delivery  port:   if  not,  the  game  is 
stopped,  and  a  "printout"  is  made  of  the  game  results  to  date  including  the 
event  that  caused  the  stoppage. 

At  one  hundred  miles  from  its  home  port,  all  ships  of  the  convoy  except 
the  forward  screen  escorts  pass  through  a  constant  barrier  (barrier  1C).   The 
losses,  if  any,  are  recorded  and  subtracted  from  the  convoy. 

After  having  transited  the  first  constant  barrier,  the  convoy  becomes 
liable  to  an  air  attack.   This  attack  is  played  as  the  probability  of  a  con- 
voy ship  being  sunk  per  day  at  sea  by  air  attack.   The  convoy  then  enters 
the  submarine  patrol  area  and  may  encounter  a  submarine  while  in  this  area. 
The  submarine/convoy  encounter  will  be  described  in  greater  detail  later. 

Ships  lost  because  of  any  one  of  the  aforementioned  reasons  are  noted 
by  type  of  ship  sunk,  method  of  sinking,  and  time  of  loss,  and  are  subtracted 
from  the  convoy  to  be  printed  out  at  the  close  of  the  simulation.   If  these 


*  Variations  in  inputs  permit  the  forward  screen  to  be  composed  alternately 
of  aircraft,  sonobuoys ,  submarines,  etc. 
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losses  include  ASW  escorts  or  forward  screen  escorts,  the  convoy's  ability 
to  detect**  enemy  submarines  will  be  diminished.   Each  time  ASW  escorts 
and  forward  screen  escorts  are  sunk,  new  lower  detection  probabilities  for 
the  various  parts  of  the  convoy  are  calculated. 

After  passing  through  the  submarine  patrol  area,  the  convoy  becomes 
liable  to  another  air  attack,  then  passes  through  a  second  constant  barrier 
one  hundred  miles  from  its  delivery  port.   After  reaching  the  delivery  port, 
the  convoy  unloads  and  stays  in  port  a  prescribed  number  of  days.   During 
this  period  the  cargo  ships  are  again  liable  to  air  attack.   Of  those  ships 
sunk  in  an  air  attack,  one-half  are  assumed  to  be  sunk  before  unloading  and 
one-half  after  unloading.   The  number  of  cargo  ships  delivered  are  now 
tabulated  and  the  convoy  is  ready  to  proceed  home. 

Proceeding  home,  the  convoy  follows  the  same  route  in  reverse  order, 
passing  through  the  same  barriers,  the  same  submarine  patrol  area,  and  the 
same  series  of  air  attacks.   When  it  reaches  its  home  port,  the  number  of 
ships  remaining  in  the  convoy  are  placed  in  a  "ship  pool"  by  ship  type. 
These  ship  pools  will  then  be  used  to  make  up  future  convoys. 

Submarines 

This  model  is  designed  to  handle  three  types  of  submarines.   All  of  the 
parameters  which  describe  the  submarine  types  are  inputs  to  the  game.   At 
the  start  of  the  war,  the  opposition  submarines,  chosen  by  type  of  submarine 
as  well  as  by  percentage  of  each  type,  are  placed  in  two  predetermined  "game 
positions,"  selected  according  to  the  option  of  the  user.   The  two  game  po- 
sitions are  "on  station"  and  "in  port."  For  those  submarines  that  start  the 
game  in  port,  a  given  number  (input)  will  leave  at  a  given  time  interval 
(input)  until  all  have  left  for  their  stations. 

While  a  submarine  is  in  port  it  may  be  susceptible  to  strikes  at  source. 
In  addition,  the  submarine  building  rate  creates  new  submarines  as  the  war 
progresses.   As  each  new  submarine  or  "turnaround"  submarine  leaves  port, 
it  is  given  a  full  load  of  anti-ASW  and  anti-shipping  torpedoes,  and  a  time 
at  which  it  must  depart  from  station  and  return  home.   Both  of  these  inputs 
are  a  function  of  submarine  type.   Submarines  that  are  on  station  at  the 
start  of  the  war  are  given  a  full  load  of  both  types  of  torpedoes,  but  their 
time  remaining  on  station  is  determined  by  a  uniform  random  distribution. 
The  enemy  submarine  patrol  areas  are  considered  to  be  a  constant  distance 
from  the  submarine  port. 

While  all  submarines  are  considered  to  transit  the  same  average  distance 
to  their  patrol  area,  they  can  appear  in  one  of  four  areas  along  the  convoy 
route,  each  area  involving  either  good  or  bad  sonar  conditions.   Each  area 
is  looked  upon  as  a  "sonar  box"  and  is  given  a  certain  percentage  of  each 
type  of  submarine  (input).   A  convoy  must  pass  through  the  sonar  boxes  in  a 


**  The  terms  "detect"  or  "detection"  as  used  throughout  this  research  contri- 
bution refer  to  the  process  of  sensing  a  signal  and  classifying  it  as  due 
to  a  target. 
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fixed  sequence,  viz.,  Box  1,  first,  Box  2,  second,  etc.   Either  good  or 
bad  sonar  conditions  may  be  assigned  to  each  box.   If  a  given  submarine  is 
put  into  a  certain  box,  the  submarine  derives  both  its  detection  and  sinking 
probabilities  with  respect  to  the  convoy,  not  only  from  its  type,  but  also 
from  the  sonar  conditions  assigned  to  the  box.   Thus  for  each  decision  point 
in  the  submarine/convoy  engagement,  there  is  in  effect  a  2  x  3  probability 
matrix,  the  rows  of  which  correspond  to  the  two  possible  sonar  conditions, 
and  the  columns  to  the  three  possible  types  of  submarines. 

When  the  submarine  leaves  port,  it  may  be  required  to  pass  through  as 
many  as  eight  barriers,  four  constant  and  four  variable,  before  arriving  in 
the  submarine  patrol  area.   Each  barrier  can  be  used  or  not  used  as  desired 
(input).   The  constant  barriers  are  of  the  same  type  as  those  described  for 
the  convoy,  except  that  in  addition  to  having  a  probability  of  being  sunk 
as  it  transits  the  barrier,  the  submarine  can  also  incur  a  time  delay  as  it 
goes  through  or  around  the  barrier.   Any  time  delay  so  incurred  will  of 
course  cut  down  the  remaining  time  the  submarine  can  spend  in  its  patrol 
area. 

When  the  submarine  crosses  a  variable  barrier,  it  not  only  has  the  chance 
of  being  sunk  itself,  but  also  has  the  chance  of  sinking  a  ship  in  the  barrier. 
If  the  submarine  does  sink  a  ship  in  the  barrier,  the  probability  of  detection 
for  this  barrier  is  assigned  a  new  lower  value.   The  barrier  ship  sunk  is 
replaced  at  some  later  time  in  the  game  (an  input),  and  the  detection  prob- 
ability of  the  barrier  increased.   In  this  way  the  variable  barrier  not  only 
destroys  some  of  the  submarines  that  must  transit  it,  but  is  itself  sub- 
jected to  possible  reductions  in  effectiveness.   Of  course,  submarines  which 
are  detected  by  the  barrier  but  not  sunk  suffer  time  delays  which  reduce 
the  time  that  can  be  spent  in  the  patrol  area. 

As  long  as  the  submarine  is  at  sea,  it  becomes  susceptible  to  being  found 
and  engaged  by  a  HUK  group.   There  are  two  different  probabilities  per  day  at 
sea  that  a  submarine  will  encounter  a  HUK  group;  one  applies  when  the  sub- 
marine is  in  transit  to  or  from  its  patrol  area,  and  the  other  when  the  sub- 
marine is  in  its  patrol  area  waiting  for  a  convoy.   Should  the  submarine  be 
found  by  a  HUK  group,  the  HUK  group  will  attempt  to  sink  the  submarine.   In 
turn,  the  submarine  has  a  chance  of  sinking  one  of  the  HUK  ships. 

The  time  of  arrival  "on  station"  depends  upon  the  speed  of  the  submarine 
(which  in  turn  depends  upon  submarine  type)  and  upon  the  time  delays  the 
submarine  undergoes  while  crossing  barriers  and  fighting  HUK  groups.   The 
submarine  remains  on  station  or  in  its  patrol  area  (sonar  box)  until  a  con- 
voy crosses  the  "subline."  The  length  of  the  "subline"  or  what  could  better 
be  called  the  submarine  patrol  area,  is  an  input. 

It  is  assumed  that  a  convoy  crosses  the  center  of  the  submarine  patrol 
line.  When  this  happens,  the  first  submarine  is  drawn  from  Box  1  and 
placed,  using  a  uniform  random  distribution,  on  the  subline.  Its  position 
on  this  line,  relative  to  the  center  of  the  line,  is  compared  to  an  athwart- 
ship  detection  range  for  this  type  of  submarine  and  for  its  sonar  box.  If 
it  is  within  the  detection  range,  the  submarine  detects  the  convoy,  other- 
wise it  does  not.  If  it  detects  the  convoy,  it  will  attempt  to  "close"  on 
this  convoy;  if  it  does  not  detect  the  convoy,  it  will  be  returned  to  its 
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sonar  box  (patrol  area). 

The  first  three  boxes  are  played  in  this  manner:   the  fourth  box  is 
played  slightly  differently.   As  the  convoy  approaches  the  fourth  box 
(fourth  patrol  area),  there  is  given  a  probability  per  month  of  the  war 
that  submarines  in  this  patrol  area  will  contact  the  convoy.   If  contact 
is  made,  this  box  will  be  played  as  the  other  three.   If,  however,  this 
contact  is  not  made,  this  box  will  not  be  played  for  this  convoy.   This 
fourth  box  might  then  be  used  for  a  submarine  wolf  pack;  for,  although  a 
wolf  pack  has  a  low  probability  of  intercepting  a  convoy,  given  this  in- 
terception it  has  a  high  probability  of  sinking  ships  in  the  convoy. 

In  general,  if  a  submarine  has  detected  a  convoy,  it  will  next  try 
to  penetrate  the  convoy  defenses  and  sink  cargo  ships.   (However,  an  option 
is  permitted,  if  desired,  by  which  submarines  will  try  to  sink  A5W  escorts 
in  an  anti-ASW  escort  campaign. )   The  event  that  plays  the  submarine  en- 
gaging the  convoy  is  Event  Subroutine  No.  19,  "Submarine  vs.  Convoy."   In 
this  event,  the  submarine  must  first  penetrate  the  forward  screen  (or  out- 
ward screen  of  the  convoy).   Given  a  successful  penetration  of  this  screen, 
it  must  then  penetrate  the  ASW  screen  (or  innermost  convoy  screen).   When 
it  has  successfully  penetrated  both  of  these  screens,  it  may  then  attack 
cargo  ships.   When  the  submarine  has  sunk  one  or  more  ships,  either  cargo 
or  screen  ships,  the  surviving  screen  ships  form  surface  attack  units 
(SAU's)  which  attempt  to  locate  the  submarine  and  attack  it.   Anti-shipping 
and  anti-ASW  torpedoes  are  used  during  this  attack,  and  whenever  the  allow- 
able minimum  number  of  torpedoes  of  the  type  needed  to  continue  the  attack 
has  been  reached,  the  submarine  is  required  to  break  off  the  attack  and 
start  its  transit  for  home.  ... 

After  a  submarine  has  either  reached  its  minimum  number  of  torpedoes 
or  has  used  its  allowable  time,  whichever  occurs  first,  it  must  return  to 
its  home  port.   On  the  way  home,  the  submarine  must  traverse,  in  reverse 
order,  the  same  barriers  that  it  crossed  on  its  way  out  to  station.   When 
the  submarine  arrives  at  its  home  port,  it  will  stay  there  for  a  certain 
number  of  days  (input)  and  then  be  resupplied  and  again  sail  for  its  patrol 
areas 
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